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demonstrate  the  feasibility  of  designing  a universal  mount  fixture, 
a three  '•’agree  of  freedom  mathematical  model  was  developed  that  simulates  the 
man  and  weapon  as  a coupled  dynamical  system.  A sensitivity  analysis  was 
performed  on  both  analog  and  digital  computers  to  determine  which  of  the 
system  parameters  are  critical  to  the  design  of  the  mount  and  to  obtain  bounds 
bn  them.  The  mathematical  simulation  was  supported  by  an  extensive  test  firing 
program,  involving  shooters  of  various  sizes  firing  the  K16  rifle  and  the 
M79  and  M203  grenade  launchers.  In  these  tests  the  shoulder  reaction  force 
and  the  dynamical  motions  of  the  man-weapon  system  were  recorded-^Inforraation 
provided  by  the  tests  was  used  to  estimate  some  of  the  critical  p^rhmeters  of 
the  system  aud  to  provide  a means  for  validating  the  model  prediction's. 

Results  of  this  study  indicate  the  feasibility  of  designing  a universal  small 
arms  mount  based  on  a two  degree  of  freedom  concept  that  would  be  capable  of 
testing  a variety  of  small  arm  weapons. 
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1.0  INTRODUCTION 


For  many  years  the  Army  has  encountered  problems  in  hard  stand  test 
firings  of  shoulder  supported  weapons.  There  are  two  basic  causes  for 
these  problems: 

a.  Test  fixtures  do  not  adequately  represent  actual  firing 
situations,  which  results  in  an  incongruity  between  data  obtained  in  test 
firings  and  actual  field  use  firings. 

b.  Test  fixtures  are  not  universal,  so  that  a particular  weapon 
might  reveal  one  set  of  characteristics  in  one  type  of  fixture  and  a 
different  set  in  a second  type  of  fixture.  For  the  above  reasons  many 
weapon  producers  have,  in  the  past,  been  forced  to  shut  down  production 
because  their  weapons  failed  to  pass  acceptance  tests.  It  is  estimated 
that  over  $100,000  could  have  been  saved  on  the  M16  rifle  program  if  only 
one  of  the  many  production  shut  downs  due  to  weapon  malfunctioning  was 
eliminated . 

The  objectives  of  this  report  are  to  determine  the  feasibility  of 
designing  a small  arms  mount,  that  adequately  simulates  the  mount  reaction 
force  for  an  actual  shoulder  supported  firing  configuration,  and  to  pro- 
pose a design  concept  of  a mount  fixture  and  arrive  at  the  critical  design 
parameters  for  the  fixture. 

To  show  the  feasibility  of  designing  a mounting  fixture,  a prototype 
mathematical  model  was  developed  that  simulates  the  man  and  weapon  as  a 
coupled  dynamical  system.  The  derivation  of  the  model  equations  is 
presented  in  Section  2.0. 

Successfully  modeling  man  as  an  integral  part  of  a weapon  system  rests 
upon  two  important  facts.  First,  the  dynamical  motions  that  occur  in  a 
typical  shooting  situation  are  characterized  by  small  angle  oscilliatory 
motions  of  the  human  operator  about  an  initial  aiming  position.  Second, 
the  applied  breech  pressure  force  consists  of  a periodic  sequence  of 
impulses.  Man’s  neuro-muscular  reaction  time  is  slow  compared  with  the 
weapon  firing  rate  and  his  ability  to  think  and  actively  respond  to  these 
force  inputs  does  not  influence  the  dynamics  of  the  system  until  after  a 
significant  amount  of  time  has  elapsed.  An  initial  passive  response  phase 
of  motion  is,  therefore,  a characteristic  of  the  man-weapon  interaction 
problem. 
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Before  listing  the  specific  modeling  assumptions,  a discussion  of  the 
exact  nature  of  human  dynamical  response  is  presented.  The  human  body  is 
a nonhomogeneous  composition  of  body  segments  with  multiple  degrees  of 
freedom.  In  addition  to  being  nonhomogeneous,  each  body  segment  is  deformable 
because  of  blood  flow  and  muscle  action;  furthermore,  the  human  neuro- 
muscular system  behaves  as  a servo  control  mechanism  to  force  inputs. 

Some  pecularities  such  as  blood  flow  have  little  affect  on  the  dynamical 
behavior  of  the  human  body.  The  dynamical  representations  of  the  body 
segments  and  the  representation  of  the  shooter’s  behavior  as  a servo  control 
mechanism  is  treated  by  making  certain  mechanical  assumptions. 

Because  of  the  overall  complexity  of  the  human  body,  mathematical 
models  of  mechanical  systems  involving  human  body  interactions  must  per- 
tain to  rather  specific  configurations.  Instead  of  modeling  the  man- 
weapon  interaction  problem  for  all  types  of  firing  positions,  a more  suitable 
analytical  approach  is  to  model  one  specific  firing  position.  In  this 
way  the  kinematical  constraints  become  easier  to  prescribe  and  the  resulting 
equations  of  motion  become  less  complex.  Usually  two  or  three  firing 
positions  are  standard  for  small  arm  weapons.  The  vertical  standing  position, 
with  the  rifle  held  up  against  the  right  shoulder,  was  chosen  for  this 
analysis  primarily  because  the  motions  of  the  man  and  weapon  in  this  position 
can  be  adequately  represented  with  only  three  independent  degrees  of  free- 
dom. 

For  convenience  the  time  history  of  motion  of  the  man-weapon  system 
is  separated  into  a set  of  distinct  phases.  In  particular,  the  human 
operator's  neuromuscular  response  is  separated  into  a passive  phase  followed 
by  an  active  phase.  Breech  pressure  forces  for  the  M16  rifle  consist  of  a 
series  of  impulses  of  approximately  one  millisecond  duration,  which  repeat 
periodically  approximately  every  80  milliseconds.  Because  the  human  body 
is  not  capable  of  reacting  actively  to  force  inputs  during  the  first 
150-200  ms,  this  initial  time  interval  is  characterized  by  a passive  response 
of  the  human  operator.  Afterwards,  the  neuromuscular  system,  which  acts 
somewhat  as  a servo  control  mechanism,  can  significantly  influence  the  wea- 
pon motion;  however  the  physical  effects  of  the  man's  active  neuromuscular 
response  generally  are  not  noticeable  until  after  300  milliseconds  have 
elapsed. 
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Mach  insight,  in  the  problem  of  modeling  man-weapon  interactions, 
was  provided  by  observations  of  high  speed  photographic  films  that 
were  obtained  from  tests  conducted  at  the  Harry  Diamond  Laboratory. 

In  these  films,  both  the  top  view  and  the  side  view  of  a man  firing  the 
M16A1  rifle  in  automatic  fire  were  recorded.  A stationary  grid  was  placed 
behind  the  test  subjects  to  provide  an  inertial  reference  system.  The  rifle 
was  held  by  the  test  subjects  in  each  of  three  separate  firing  positions: 

a.  off  the  hip 

b.  on  the  right  shoulder  with  the  body  in  a vertical  standing  position 

c.  on  the  right  shoulder  with  the  body  in  a slanted  position 

The  vertical  standing  position  was  selected  for  the  theoretical  analysis. 

Observations  of  the  films  helped  to  establish  the  predominant  degrees 
of  freedom  of  the  system  and  the  time  history  of  motion  in  each  coordinate. 
Additional  high  speed  photographic  data  were  obtained  by  Mr.  Thomas  Hutchings 
for  the  M79  and  M203  grenade  launchers.  This  data  provided  the  time  history 
motions  for  the  higher  impulse  weapons  (the  M79  and  M203  grenade  launchers 
generate  an  impulse  of  about  2.5  lb-sec  compared  with  1.2  lb-sec  for  the 
M16  rifle).  Additional  data  was  obtained  for  the  M16  rifle  at  the  Keith 
L.  Ware  Simulation  Laboratory,  Rock  Island  Arsenal.  A load  cell  device 
was  used  to  measure  the  transmitted  shoulder  mount  force  and  two  displacement 
transducers  were  used  to  obtain  the  weapon  rotation.  These  experiments 
are  discussed  more  fully  in  Section  3.0.  The  data  from  these  experiments 
provided  considerable  insight  on  how  to  construct  the  analytical  model 
and  also  provided  a means  for  estimating  some  of  the  unknown  system  parameters. 

Section  4.0  contains  the  results  of  a sensitivity  analysis  that  was 
performed  with  the  man-weapon  interaction  model.  As  a result  of  this 
analysis,  two  prototype  small  arms  mount  designs  are  recommended.  The 
conclusions  and  recommendations  of  this  study  are  presented  in  Section  5.0. 
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2.0  DERIVATION  OF  THE  BIOMECHANICAL  EQUATIONS  OF  MOTION  FOR  SHOULDER 

FIRED  SMALL  ARMS 

2.1  Test  Fixture  Design  Procedure 

The  design  specifications  for  the  universal  small  arms  mount  fixture 
were  obtained  from  analytical  simulations  of  the  biomechanical  inter- 
1 action  problem.  That  is,  the  design  configuration  and  design  parameters 
of  the  mount  are  based  on  results  of  a dynamic  analysis  of  the  man-weapon 
interaction  problem.  The  specific  steps  followed  in  this  procedure  are 
listed  below: 

a.  Determine,  through  observations  of  high  speed  film  data,  the 
minimum  required  number  of  independent  degrees  of  freedom  needed  to 
represent  the  biomechanical  interaction  forces  for  shoulder  supported 
small  arm  weapons. 

b.  Develop  a biodynamical  model  of  the  man-weapon  system  and 
determine  the  sensitivities,  of  mount  force  and  of  the  pitching  motion  of 
the  weapon,  to  variations  in  the  various  biomechanical  system  parameters. 

c.  Obtain  a preliminary  design  for  a universal  small  arms  test 
fixture  based  on  the  biomechanical  analog  and  obtain  bounds  on  the 
critical  design  parameters  from  results  of  the  sensitivity  analysis. 

2 . 2 High  Speed  Photographic  Test  Results 

In  order  to  determine  the  number  of  independent  degrees  of  freedom 
that  are  associated  with  the  man-weapon  interaction  problem,  several  high 
speed  photographic  film  tests  were  performed  for  shoulder  supported  firings 
of  the  M16A1  rifle  and  the  M79  and  M203  grenade  launchers.  These  tests 
Involved  several  shooters  of  various  weights  and  heights  so  that  the  effects 
of  variations  in  human  mount  characteristics  would  be  observed.  The  data 
extracted  from  these  tests  include  the  transmitted  recoil  force  and  motions 
of  the  weapon  and  the  shooter.  The  shoulder  support  firing  configuration 
was  selected  for  the  tests,  because  the  weapon  mount  forces  are  easier  to 
measure  and  the  biodynamical  motions  are  less  complex  than  for  other  standard 
firing  configurations,  such  as  the  off-hip  configuration. 
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Observations  of  the  high  speed  photographic  film  data  led  to  the 
development:  of  a man-weapon  biodynamics  model.  The  films  recorded  the 
test  subjects  in  both  the  top  and  side  views.  The  top  view  was  obtained 
by  a 45°  inclined  mirror  located  directly  above  the  test  subject.  During 
the  tests,  the  shooters  fired  five  rounds  from  the  M16A1  rifle  in  single 
shot,  semi-automatic,  and  burst  modes,  and  fired  single  shots  from  the 
grenade  launchers  at  various  range  settings.  Observations  of  the  film 
sequences  for  these  tests  indicated  which  rotational  and  translational 
motions  of  the  weapon  and  shooter  are  predominant. 

In  addition  to  the  translatory  motion  of  the  weapon  toward  the 
shoulder,  two  predominant  motions  of  the  man-weapon  system  were  observed. 

These  motions  include: 

a.  A rotation  in  the  vertical  plane  of  the  upper  torso  of  the  man, 
who  initially  (for  approximately  0.3  second  duration)  pivoted  about  his  hips. 
Afterwards  the  man  appears  to  consciously  react  to  the  weapon  recoil  and 

his  motion  becomes  more  complex. 

b.  A rotation  in  the  vertical  plane  of  the  weapon  and  the  shooter's 
arms,  pivoting  about  the  shoulder.  The  top  view  revealed  a negligible 
amount  of  yawing  of  either  the  shooter  or  the  weapon  compared  with  the 
vertical  pitching  rotations.  Selected  angles  and  position  coordinates 
were  measured  for  each  film  sequence  using  the  Vanguard  Motion  Analyzer 

at  the  University  of  Iowa  Hospital  Biomechanics  Laboratory  in  Iowa  City,  Iowa. 

Results  of  the  high  speed  photographic  film  tests  reveal  a period 
of  relative  inactivity  in  the  pitching  rotations  for  approximately  the 
first  20-50  milliseconds  after  the  commencement  of  each  shot.  This  delay 
effect  is  probably  caused  by  the  compression  of  weapon  padding  and  soft 
body  tissue  as  the  weapon  translates  toward  the  shoulder.  After  the  weapon 
is  fully  compressed  against  the  shoulder  the  remaining  kinetic  energy  of 
the  weapon  induces  the  rotational  pitching  motions  of  both  the  man  and  the 
weapon.  Yawing  motion  apparently  was  minimized  by  the  stance  of  the  man, 
as  the  weapon  is  held  almost  parallel  with  the  breadth  of  his  chest  during 
firing.  The  net  torque  produced  about  the  man's  vertical  yaw  axis  is 
therefore  small. 
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2.3  Man-Weapon  Model  Assumptions 

A mathematical  model  simulating  the  man-weapon  biodynamical  interactions 
from  weapon  recoil  was  then  developed  based  on  information  supplied  by  the 
film  data  and  on  biomechanical  data  supplied  by  various  sources.  Figure  1 
contains  a schematic  representation  of  the  model.  The  X^,  X^,  X^ 

coordinate  system  is  a fixed  system  with  the  origin  at  the  man's  hip. 
Coordinate  system  Y^,  Y^,  Y^  is  attached  to  the  right  shoulder  pivot  and 
rotates  with  the  weapon.  The  Y^  axis  is  parallel  to  the  gun  barrel  center- 
line.  Variable  x locates  the  center  of  mass  of  the  weapon  combined  with  the 
man's  arms.  Variable  0 measures  the  absolute  pitch  of  the  man  in  the  - 
X^  plane  and  variable  <|>  measures  the  pitch  of  the  weapon  relative  to  the 
man’s  trunk.  The  total  pitch  of  the  weapon  is  therefore  the  sum  of  angles 
0 and  <j>.  Xo 


Figure  1 Man-Weapon  Interaction  Model 
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Because  the  development  of  the  biodynamics  involves,  in  part,  a 

mechanical  description  of  the  human  body,  several  simplifying  assumptions 
were  required  in  order  to  arrive  at  a fairly  simplified  but  realistic 
representation  of  the  system.  There  is  no  simple  or  easy  approach  to 
modeling  human  body  dynamics  and  very  little  useful  data  is  available  on 
biomechanical  properties  of  the  human  body.  The  feasibility  of  successfully 
modeling  man  as  an  integral  part  of  a weapon  system  rests  upon  several 
important  considerations,  which  are  discussed  below: 

(1)  The  breech  pressure  force  for  small  arm  weapons  - rifles  and 
grenade  launchers  - is  impulsive.  A typical  breach  pressure  force  curve 
for  the  M16A1  rifle  is  illustrated  in  Figure  2 . Normally  the  M16A1 
generates  approximately  1.2  lb-sec  impulses  of  about  one  millisecond 
duration  that  repeat  periodically  every  80  milliseconds. 

(2)  During  the  initial  portion  of  a burst,  the  shooter  responds 
passively  to  transmitted  mount  forces.  Thus,  even  though  he  may  anticipate 
the  impulsive  recoil  force  and  prepares,  his  muscles  accordingly,  there  is 
an  initial  time  interval  during  which  the  shooter  is  unable  to  respond 
actively  to  the  pulse. 

The  duration  of  this  passive  phase  depends  primarily  on  the  shooter's 
neuromuscular  reaction  characteristics.  Moreover,  the  actiye  response 
of  the  shooter  intuitively  should  not  have  much  influence  on  the  time 
history  of  the  transmitted  recoil  force.  One  might  expect,  however, 
that  the  active  human  control  response,  after  approximately  0.3  seconds 
or  so,  might  have  a significant  effect  on  the  angular  positioning  of  both 
the  weapon  and  the  man's  upper  torso. 

c.  Observations  of  the  high  sp  ad  photographic  film  data  reveal 
that,  during  the  passive  response  phase,  the  shooter's  hips  remain 
stationary;  consequently  only  the  upper  torso  and  arm  segments  are 
represented  in  the  biodynamics  model. 

d.  High  speed  photographic  films  reveal  that  the  maximum  angular 
pitching  motions  of  the  man's  upper  torso  and  the  weapon  are  generally  less 
than  10°.  A small  angle  approximation  is,  therefore,  applicable  in  the 
formulation  of  the  biodynamical  equations  of  motion.  Small  angle  approximations 
are  particularly  useful  for  generating  problem  solutions  on  the  analog 
computer . 
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Figure  2 M16  Breech  Force  vs  Time 
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e.  The  human  body  can  be  represented  dynamically  as  a systei  >f  linked 
rigid  bodies  that  are  internally  stable  and  homogeneous;  furthermore,  each 
rigid  body  can  be  represented  by  simple  geometric  forms  (i.e,  elliptical 
cylinder,  frustrum  of  a cone,  etc).  These  approximations  of  the  human  body 
have  been  used  in  the  works  published  by  the  investigators  listed  in 
References  1 3 . 

2.4  Historical  Background  on  Human  Biodynamlcal  Research 

The  lumped  mass  approach  to  modeling  the  dynamics  of  the  human  body 
requires  the  specification  of  segment  masses,  moments  of  inertia  and  centers 
of  mass.  Active  interest  in  determining  the  segment  characteristics  of  the 


*Wlli  tootl  , CM'..,  "Soiii.'  Dyitu,..!.-  R>'S|ion»r  Chill  ilrtcr  inf  i nil  of  Wei  J*.lit.li'HH 
Man,"  MuhUm  of  {.cl once  ’J In  •••Ik , All  force  Institute  «f  Technology, 

Wt  I I'.liL- I’aM  <■■•).  mm  Air  force  i'.i  rrti , Ohio,  AM''!,-  TK-o'l-] 0 , All  A J 2 'j  A J , 1-962 

2McCr;ii,k,  J.M.  <md  Soger,  "Torque  Frw:  Kol  ;il  iona I Dynamics  of 

n V.uUihle  Couf  i.;,,iu  at  Lou  body  (Application  l.o  WoJghtlone  Man),"  Hunter 
of  Science  Tlirsii.,  Air  force  inr  finite  of  Technology,  Wright-fat rernon 
Air  I'm  ce  rare,  Ohio,  AD  0lO!’.l9,  19  6 A 

3Mcllenry,  11. R.  and  Knob,  K.N.,  "Coiiuuiter  Simulation  of  Urn  Automobile 
Urnuh  Victim  in  a frontal  Col  Union — A Validation  Study,"  Cornell  Aeronauti- 
cal laboratory , Inc.,  lluffolo,  Hew  York,  CAl.  Report  Mo.  4B-212G-V-1R,  I960 
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human  body  has  been  undertaken  by  numerous  investigators.  Braune  and  Fisher 
in  1889  performed  experimental  research  on  three  cadavers  and  published 
a comprehensive  study  of  weights,  volumes,  and  centers  of  mass  of  the  body 
and  its  segments.  Later,  in  1906,  Fisher  5 experimentally  determined  the 
moments  of  inertia  of  segments  from  a single  cadaver.  Since  that  time 
interest  in  this  subject  diminished  until  rather  recently.  Dempster  6,  in 
1955,  published  a study  on  human  biomechanics  that  contains  experimental  data 
taken  from  eight  cadavers. 

The  specific  data  presented  in  this  study  includes  segment  weights, 
moments  of  inertia,  densities,  center  of  mass  locations  and  volumes.  Using 
data  compiled  by  Dempster,  Braune  and  Fisher,  Barter7  in  1957  prepared  a 
series  of  regression  equations  for  predicting  body  segment  weights  from  total 
body  weight.  These  equations  have  been  used  extensively  by  biomechanical 
engineers  and  designers.  In  1963  Santsehi0  and  his  co-workers  reported  on 
moments  of  inertia  and  centers  of  mass  of  sixty-six  live  test  subjects. 
Santsehi  attempted  to  answer  the  pertinent  question  of  whether  or  not  body 
segment  parameters  can  be  predicted,  to  a reasonable  degree  of  accuracy, 
from  anthropometric  dimensions.  He  discovered  a high  correlation  factor 
of  segment  centers  of  mass  and  moments  of  inertia  with  an  individual's 
anthropometric  dimensions. 

2.5  System  Kinetic  Energy 

The  formulas  that  were  used  to  predict  the  physical  characteristics 
of  the  man's  upper  torso  and  arm  segments  are  listed  in  Appendix  A.  Since 
the  man-weapon  interaction  problem  has  been  reduced  to  three  independent 
degrees  of  freedom,  the  segment  masses  corresponding  to  the  man's  upper 
torso  (i.e.  the  head  and  torso)  are  combined  to  form  a single  rigid  body 


‘‘Braune,  W.  mid  Fincher,  0.,  "The  Center  of  Gravity  of  the  Human  Body 
no  Belated  to  the  German  Infantryman Leipzig,  ATX  138432,  1809 

®l'j  Rchcr , 0.  , "Theoretic..!  Fundament  a In  for  a Mechanics  of  hiving  Bodies 
with  8pcfci.nl  Appl icnl.JoiiH  to  Man  as  Well  an  to  Somi  Processes  oT  Motion 
of  Machines,"  11.(1.  Tubuer,  BerJ  in , ATI  153668,  1900 

^Dempster,  W.T.,  "lipnee  Requirements  of  the  Seated  Operator,"  Wright: 

Air  Development  Center,  TU-35-159,  Wright-I'at terson  Air  Force  Base, 
Ohio,  AD  87892,  1955 

7Uartor,  J.T.,  "Hst  i ration  of  the  Mass  of  tody  Scgi..entn,"  Wright  Air 
Development  Center,  Wright-Patterson  Air  Force  llano,  Ohio,  TR-57-260, 

AD  118222,  1957  ' 

°Snntsehi,  W.K.,  Dulloio,  J.,  and  Omoto,  C.  , "Moments  of  Inertia  mid 
Centers  of  Gravity  of  the  Living  Human  Body,  Aerospace  Me.dieal 
Research  tabornl  ories,  Wright-  ratterson  Air  Force  Base,  Ohio, 

AMRt -TDK-63-00,  AD  410451,  1963 
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and  the  segment  masses  corresponding  to  the  man's  arm  segments  and  hands 
and  the  weapon  are  combined  to  form  a second  rigid  body.  The  rigid  body 
comprising  the  trunk  and  the  head  is  allowed  to  rotate  in  the  vertical 
plane  about  the  man's  hip,  while  the  rigid  body  comprising  the  man's  arms 
and  the  weapon  is  allowed  to  translate  in  the  direction  of  the  recoil  force 
and  rotate  vertically  about  the  shoulder  hinge.  Any  bending  of  the  man's 
aims  during  a firing  event  is  considered  small  and  is  neglected. 

Referring  to  Figure  1,  the  kinetic  energy  of  the  system  is  given  by 


* [ 


• 2 

I 0+1 
cm  r 


(6  + <p)2  + li  r 


cm 


t + m r 
cm  r 


(2-1) 


where 

I * Moment  of  inertia  about  the  center  of  mass  of  the  head-trunk 
CU  rigid  body. 

Ir  * Moment  of  inertia  about  the  center  of  mass  of  the  arm-weapon 
rigid  body. 

M ■ Mass  of  the  head-trunk  rigid  body. 

mr  ■ Mass  of  the  arm-weapon  rigid  body. 


r 


cm 


Time  derivative  of  the  center  of  mass  location  of  the  head-trunk 
rigid  body,  and, 


r =»  Time  derivative  of  the  center  of  mass  location  of  the  arm-weapon 
rigid  body. 

In  the  X^,  Xj*  frame,  the  components  of  vectors  rcm  and  r are 


rcm  (~sin  0 + cos  0 *3)  (2-2) 

[x  cos  (0  + <f>)  - L sin  0]  X2 

[x  sin  (0  + <j>)  + L cos  0]  xj 

Taking  the  time  derivatives  of  expressions  (2-2) 
into  equation  (2-1)  leads  finally  to  result 

T " {t  (Xcm  + h t2cm)  'QZ  + Xr  + ^ 

+ m^  [ (x  - h 0 cos  <j>)2  + (x  <j>  + (x  + L sin  <f>)  $)2  (2-4) 


(2-3) 

and  (2-3)  and  substituting 
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9.6  System  Potential  Energy,  Dissipative  Mechanisms,  And  Generalized  Forces 
In  order  to  obtain  an  expression  for  the  system  potential  energy, 
the  passive  response  of  the  man's  muscles  must  be  taken  into  account. 

Prior  to  firing,  the  shooter’s  muscles  are  in  a state  of  preloading  to 
support  the  weapon  and  to  balance  the  forces  of  gravity  on  both  the  man  and 
weapon.  After  the  commencement  of  firing,  the  approximation  is  made  that 
passive  restoring  forces,  produced  by  the  various  muscles  in  the  upper  torso 
and  arms,  can  be  represented  by  a system  of  linear  springs  and  dashpots. 

This  approximation  is  justified  based  on  the  small  orders  of  magnitude  of 
the  angular  deformations  resulting  from  weapon  recoil  torques.  Torsional 
springs  and  dashpots  are,  therefore,  included  in  the  model  to  resist  motions 
in  each  of  the  two  rotational  degrees  of  freedom.  The  rearward  translatory 
motion  of  the  weapon  initially  causes  the  fleshy  padding  on  the  shoulder  to 
compress.  Any  further  rearward  displacement  of  the  rifle  results  in  the 
stretching  of  the  muscles  and  ligaments  that  attach  to  the  bone  structure 
of  the  shoulder.  Resistance  offered  by  these  muscles  and  ligaments  is 
modeled  by  a relatively  strong  spring  and  dashpot,  while  the  initial  soft 
compression  of  the  shoulder  is  modeled  by  a relatively  weak  spring.  Having 
represented  the  elastic  and  viscous  damping  characteristics  of  the  model, 
the  system  potential  energy,  V,  then  becomes 


1 k 

2 0 


(0  - 0e)2  +•§  4 - 4>e)2 
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(X  - xe)' 


(x  - xs)2  + kx 
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for  x > x„ 
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S - xe>  tx  "2  (xs  + xe)]  for  x < Xs 
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cm 


(cos  0 - cos  0q)  + m g [L(cos  0 - cos  Qq) 


+ x sin  (0  + <j>)  - Xo  sin  (Oq  + <(>q)] 


(2-5) 


where 

kg,  k,,  k^  , kx  " system  spring  constants 
w s> 

0o,  $ , xQ  « initial  conditions  for  variables  0,  <j>,  and  x respectively 
0 <f>e>  xp  = static  equilibruim  spring  positions 

g = gravity  constant 
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and 


x ■ distance  from  the  shoulder  at  which  the  shoulder  spring  stiffness 
s 

changes  value. 

Because  the  only  dissipative  mechanisms  in  the  model  are  viscous  dampers, 
a dissipative  force  potential,  P,  can  be  defined  and  is  given  by 

p ' \ + I ce  52  4 i *2 

where  cx>  cQ,  and  c^  are  the  damping  constants  for  the  variables  x,  0#  <j>, 
respectively. 

Referring  to  Figure  3,  the  generalized  forces  produced  by  the  breech 
pressure  force,  F(t),  are 

Qx  - - F(t)  (2-7*) 

— A 

Qe  “ <?D/A  x f)  * X1 

and 

% - <5d/b  * • h 

where 

rD/A  - [-L  sin  0 + x cos  (0  + <j>)  - 5 sin  (0  + <(>)]  x2 

4-  [L  cos  0 + x sin  (0  + $)  +6  cos  (0  + <j>)  ] X3 

?d/b  " x y*  + 5 ^ 

F - -F(t)  y2 


-F(t)  [cos  (0  + 4> ) x2  + sin  (0  J-  <j> ) X3I 


Expanding  the 

above  expressions  for  Q.  and  0, 

u v9 

yields  finally  the  results 

and 

% 

*■  (<$  + L cos  <£)  F(t) 

(2-7b) 

% 

- <5‘F(t) 

(2-7 c) 

r^ 

cm 

Dynamical 

Equations  of  Motion 

lent 

The  dynamical  equations  of  motion  of  the  man-weapon  interact ic.u  pro’ 
were  derived  using  the  Lugrangrl an  formulation. 
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where  L = T - V is  the  Lagrangian  potential  function.  Substituting 
equations  (2-1)  to  (2-7)  into  equations  (2-8)  and  rearranging  terms  leads 
finally  to  the  following  matrix  equation: 


where 


aH  ”a2 1 0 

~a2 1 a22  + a3  3 a23 

0 a23  a33 


in  L cos  d> 
r T 


1^,  + Mr2  + tn  L2  + 2m  Lx  sin  <j> 
cm  cm  r r 


a3  3 + m.f  L x sin  <j> 


lr  + "r  x' 


F(t)  -f  tnr  x (0  + 40  2 + m 02  L sin  <j> 


Cv.  x - ra  g sin  (0  + w) 


(2-9) 


kx  (x  “ xe>  > 
w 

kx  (x  - x8)  -I-  (xs  - xe)  , 
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for  x > x 


for  x < y. 


fg  “ <6  + L cos  <p]  F (t)  - c.Q  6 

- 2 mr  x (0  + |)  (x  + L sin  $) 

■ mr  L x J (20  + <j>)  cos  $ - kg  (0  - 6 ) 

+ M g r sin  t)  + mr  g [L  sin  0 - x cos  (0  + ifi)] 


computer  by  the  Runge-Kutta  method.  in  order  to  apply  the  Runge-Kutta 
method,  equation  (2-9)  is  premultiplied  by  the  inverse  of  the 
coefficient  matrix  and  the  resulting  system  of  three  second  order  differential 
equations  is  transformed  into  a system  of  six  first  order  differential 
equations.  This  results  in  the  set  of  matrix  equations  given  below: 
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A =>  det  [a] 
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and 


h33  “ ar  [Icm  + M rcm2  + Xr  + mr  (x  + L sin  +)2 1 * 

The  solutions  to  equations  (2-10)  were  obtained  by  separating  the  weapon 
cycle  into  two  time  intervals.  For  the  one  millisecond  time  interval,  on 
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which  the  impulsive  breech  pressure  force  is  applied,  equations  (2-10)  are 
integrated  with  the  assumption  that  the  positional  coordinates  are  con- 
stant. Thus  let 
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where  t^=  starting  time  for  the  shot  and  t^  + ^ + 0.001.  The 

integration  of  equations  (2-10)  yields  a set  of  initial  conditions  for  the 


start  of  the  second  time  interval. 
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For  the  remainder  of  the  weapon  cycle  (i.e.  until  the  commencement  of  the 
next  shot),  the  equations  (2-10)  are  solved  by  the  Runge-Kutta  method  (with 
F(t)  = 0)  and  equations  (2-11)  provide  the  initial  conditions.  This  sol- 
ution procedure  is  then  repeated  for  each  of  the  remaining  shots  in  the 
burst.  A listing  of  the  computer  program  used  to  obtain  the  numerical 
solutions  to  the  problem  is  provided  in  Appendix  B. 

ihe  digital  computer  solutions  co  the  biodynamics  problem  were  obtained 
for  a range  of  biomechanical  parameters  co ’-responding  to  a small,  average 
and  large  shooter,  and  for  a range  of  weapon  parameters  corresponding  to 
the  M16A1  rifle  and  the  M79  and  M203  grenade  launchers.  A sensitivity 
analysis  was  performed  to  obtain  estimates  of  the  unknown  spring  and  dash- 
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pot  parameters,  for  several  of  the  weapon  and  shooter  configurations.  The 
results  of  this  analysis  are  discussed  in  Section  4.0. 

In  the  second  solution  method,  the  system  of  equations  (2-9)  were 
first  simplified  by  small  angle  approximations,  scaled  and  then  programmed 
on  an  analog  computer.  Analog  techniques  are  suitable  for  performing  a 
parametric  sensitivity  analysis,  because  each  system  parameter  can  be 
conveniently  changed  by  adjusting  a pot  setting  and  the  solutions  are 
displayed  immediately  on  a recording  device.  Further  details  regarding  the 
analog  simulation  are  provided  in  Appendix  C. 
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3.0  MAN-WEAPON  INTERACTION  EXPERIMENTS 


3.1  Description  of  Experimental  Teats 

In  order  to  gather  the  critical  information  needed  to  analytically 
describe  the  man-weapon  interaction  problem,  three  separate  series  of 
experiments  were  performed  that  involved  test  firings  with  several  shooters 
of  various  sizes  and  various  weapon  types.  Data  gathered  ftom  these  tests 
include : 

a*  M16  data  obtained  from  the  Harry  Diamond  Laboratory  using  high 
speed  photographic  equipment  to  record  time  histories  of  rotations  and 
displacements  of  body  segments  and  of  the  weapon. 

b.  M79  and  M203  grenade  launcher  data  obtained  at  the  Harry  Diamond 
Laboratory  test  facility  using  high  speed  photographic  equipment  to  record 
rotations  and  displacements  of  the  weapon  and  the  man. 

c.  M16  data  from  the  Keith  L.  Ware  Simulation  Laboratory,  obtained 
using  a fixture  between  the  shoulder  and  the  weapon  to  measure  force 
reactions  at  the  shoulder  and  two  displacement  transducers  located  at  the 
front  and  rear  of  the  weapon  to  measure  rotations  of  the  weapon. 

In  both  photographic  data  experiments  at  the  Harry  Diamond  Laboratory, 
the  shooters  stood  in  front  of  a background  of  vertical  and  horizontal 
grid  lines.  These  grid  lines  formed  a fixed  reference  frame  from  which 
relative  motions  of  the  man  and  weapon  were  measured.  Both  top  and  side 
views  were  photographed.  The  top  view  was  recorded  by  use  of  an  inclined 
mirror  placed  directly  above  the  test  subject  (see  Figure  4).  The  film 
speed  of  1100  frames/second  was  sufficiently  fast  to  record  the  ejection 
of  muzzle  smoke  after  each  shot  from  the  M16  rifle.  Thus  the  time  of 
occurrence  of  each  shot  was  measured  for  the  five  shot  burst.  Fixed  reference 
points  were  selected  at  the  intersections  of  stationary  grid  lines  on  the 
side  and  top  views.  For  each  frame,  the  positions  of  several  selected 
points  on  the  man  and  weapon  were  measured  with  respect  to  the  selected 
reference  points  using  a Vanguard  Motion  Analyzer.  A computer  pro- 
gram was  then  used  to  process  and  convert  the  raw  data  into  values  of 
angular  displacements  of  the  weapon  and  man. 
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3. 2 M16A1  Test  Results 

Some  of  the  important  observations  and  interpretations  of  the  M16 
data  are  discussed  below.  Two  major  angular  motions  of  the  weapon  system 
were  observed  in  addition  to  the  translatory  motion  of  the  rifle  towards 
the  shoulder.  They  are: 

a.  An  angular  deflection  of  the  man's  trunk  that  pivoted  about 
his  hip.  This  mo Lion  is  referred  to  as  pitch  of  the  man. 

b.  An  angular  pivoting  of  the  arms  and  rifle  about  the  man's  should': r. 
This  motion  is  referred  to  as  pitch  of  the  rifle. 

A plot  of  the  two  pitch  coordinates  versus  time  is  shown  in  Figure?. 

A period  of  relative  inactivity  is  noticable  in  the  pitch  of  the  man  Just, 
after  the  first  and  second  shots.  During  this  period  of  inactivity,  the 
weapon  most  likely  is  compressing  against  the  soft  muscular  portion  of:  the 
shoulder.  A period  of  activity  in  the  pitch  coordinates  follows  the 
completion  of  the  soft  compressive  phase.  Yawing  of  the  man  and  rifle  was 
small  compared  with  the  pitching  rotations.  Yaw  was  minimized  by  the 
rigid  stance  of  the  man  and  his  initial  positioning  of  the  rifle. 

3. 3 M79  and  M203  Grenade  Launcher  Tests 

The  grenude  launcher  tests  were  performed  on  an  experimental  setup 
similar  to  the  M16  experiment  (see  Figure  6)  • One  difference  was  that 
fixed  points  on  the  man  and  the  weapon  were  identified  by  circular  pieces 
of  reflective  tape.  The  tape  markings  helped  to  simplify  the  data  reduction 
on  the  Vanguard  Motion  Analyzer.  Another  modification  in  the  grenade 
launcher  tests  was  that  the  roturary  prism  type  camera  was  replaced  with  a 
shutter  typo  camera,  operating  at  a slower  300  frames/second,  to  improve 
the  film  resolution.  The  reduced  film  speed  is  sufficient  to  ident'fy  when 
the  grenade  exits  the  barrel  and  to  obtain  continuous  tire  history  data  of 
all  significant  motions.  Several  film  sequences  were  recorded  with  different 
size  men  who  fired  from  both  the  shoulder  and  hip  configurations.  The  data 
were  reduced,  on  the  Vanguard  Motion  Analyzer,  by  the  Intertech  Corporation 
in  Iowa  City,  Iowa.  In  reduced  form,  the  data  provided  pitch  and  yaw 
rotations  of  the  weapon  along  with  the  pitch  of  the  man  and  the  weapon 
displacement . 


3-3 


Experimental  Data  For  The 


The  M203  grenade  launcher  data  were  the  most  thoroughly  studied  data. 

The  M2Q3  grenade  launcher  is  mounted  below  the  barrel  on  the  M16  rifle, 
thus  causing  this  weapon  to  pitch  downward  rather  than  upward,  as  most 
other  weapons  do.  Also  the  data  for  the  shoulder  support  configuration 

reveal  that  the  pitch  rotation  of  the  weapon  varies  somewhat  with  the 
size  of  the  shooter.  There  is  less  pitch  of  the  weapon  for  a heavier  man 
than  for  a lighter  man.  This  is  expected  because  the  heavier  man  has  greater 
inertia.  Firings  from  the  hip  produced  less  weapon  pitch  than  firings  from 
the  shoulder.  Again,  this  is  expected  since  the  constraints  on  the  weapon 
in  the  hip  firing  position  tend  to  limit  the  amount  of  weapon  pitch  rotation. 

Plots  of  the  total  weapon  rotation  versus  time  have  been  generated  for 
the  test  firings  of  the  M79  and  M203  grenade  launchers.  Several  of  the  plots 
for  the  M203  grenade  launcher  tests  appear  in  Figures  7 to  13.  The  test 
conditions  for  these  film  sequences  are  listed  in  Table  1.  Each  plot  con- 
tains both  the  measured  weapon  pitch  angle  data  points,  connected  by  straight 
line  segments,  and  a fourth  order  least  square  polynomial  approximation  to 
the  data.  The  time  interval  begins  at  the  commencement  of  firing  and  ends 
approximately  when  the  man's  active  neuromuscular  response  starts  to  influence 
the  motions  of  the  system  (i.e.  at  roughly  0.3  second).  In  each  film  sequence, 
the  weapon  pitches  downward  from  the  initial  aim  elevation  angle  to  a maximum 
value  of  about  seven  degrees.  The  shapes  of  the  plots  vary  somewhat  both  for 
different  shooters  and  different  aim  elevations.  The  plots  in  Figures  11  and 
12  reveal  that  different  test  results  were  obtained  even  for  teste  repeated 
under  identical  condi tions.  This  may  be  due  partly  because  of  errors  in 
data  reduction  and  because  the  shooter  does  not  prepare  himself  in  exactly 
the  same  way  for  eacli  shot. 


TABLE  1 M203  GRENADE  LAUNCHER  TEST  CONDITIONS 


FILM 

SEQUENCE 

SHOOTER 
WEIGHT  (LB) 

ELEVATION 
ANGLE  (DKG) 

6 

?. 

HO  . 

5* 

0 

4 

HO 

10* 

U 

6 

140 

20° 

9 

3 

103 

5* 

9 

5 

103 

:io° 

9 

6 

165 

10* 

7 

165 

20°  - 2!i° 

- 
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Figure  11  Film#  9,  Seq,#  5 - Weapon  Rotali 


Weapon  Rotation 


3.4  Measurement  of  the  Shoulder  Reaction  Force 

Additional  M16  data  were  obtained  at  the  Keith  L.  Ware  Simulation 
Laboratory.  The  experimental  setup  is  shown  in  Figure  14.  A shoulder 
fixture  (see  Figure  15)  containing  load  cells  and  an  accelerometer  was 
used  to  measure  the  weapon  recoil  force  transmitted  to  the  man's  shoulder. 

The  force  reaction  at  the  shoulder  is  obtained  by  subtracting  the  inertial 
force  of  the  fixture  from  the  measured  load  cell  force  (see  Figure  16). 

The  linear  potentiometers  measured  the  vertical  displacement  of  the  weapon 
at  two  points  on  the  rifle.  Referring  to  Figure  17,  the  weapon  pitch 
rotation,  8,  is  then  calculated  from  the  difference  between  the  two  vertical 
displacements . 

The  most  significant  conclusion  that  was  drawn  from  the  Ware  Laboratory 
data  is  that  the  shoulder  reaction  force  is  of  very  short  duration  (i.e. 
impulsive)  compared  with  the  rifle  firing  cycle  time.  The  maximum  forces 
ranged  from  20  to  180  pounds.  Peak  recoil  forces  tended  to  be  greater  for 
the  larger  riflemen  than  for  the  smaller  riflemen.  Also,  the  peak  values 
tended  to  decrease  significantly  from  round  to  round.  The  ratios  of  peak 
values  for  successive  shots  varied  considerably  among  the  tests.  Ratios 
for  successive  peak  recoil  forces  in  the  neighborhood  of  0.85,  however,  were 
common. 

The  data  gathered  from  these  experiments  provided  insight  for  constructing 
the  simulation  model  and  provided  a means  for  estimating  some  of  the  unknown 
parameters  in  the  model  (i.e.  stiffness  and  damping  constants).  The  time 
history  data  obtained  by  high  speed  photography  is  more  reliable  than  the 
data  measured  by  linear  potentiometers,  because  the  photographic  technique 
does  not  have  any  direct  influence  on  the  experiment.  On  the  other  hand, 
the  use  of  high  speed  photography  to  measure  motions  has  the  disadvantage  of 
requiring  a large  amount  of  time  and  effort  in  data  reduction. 


3-14 


SHOULDER  FORCE  MEASURING  DEVICE 

My  - MASS  OF  FIXTURE 

~ MEASURED  FORCE  (LOAD  WITH  CELLS) 

X - MEASURED  ACCELERATION 

rs  - COMMUTED  FORCE  EXERTED  ON  THE  SHOULDER 

h “ *rL  ” *** 


Figure  16  Computation  of  Transmitted  Recoil  Force 


BIODYNAMIC AL  SENSITIVITY  ANALYSIS 


4.0 

4.1  Introduction 

A sensitivity  analysis  was  performed  to  determine  which  biomechanical 
parameters  of  the  man-weapon  system  are  critical  to  the  design  of  a universal 
small  arms  mount.  In  this  analysis,  both  the  peak  force  at  the  shoulder 
and  the  absolute  weapon  pitch  rotation  after  0.3  seconds  were  taken  as  the 
primary  output  variables.  The  sensitivities  of  these  output  variables  to 
systematic  variations  in  each  of  the  system  parameters  were  obtained  by 
numerical  and  analog  computer  solutions  of  the  biedynamical  equations  of 
motion. 

Both  single  and  burst  fire  modes  were  analyzed.  The  driving  force 
F(t)  in  equations  2-7  represents  what  is  commonly  called  the  breech  force. 

The  important  feature  of  this  function  is  the  short  duration  (about  one 
millisecond)  compared  with  the  total  cycle  time  (approximately  80  milli- 
seconds) . For  the  digital  computer  simulation,  the  function  F(t)  is 
approximated  by  Dirac's  delta  function  — a symbolic  function  commonly  used 
in  mathematical  physics  that  has  finite  area  but  is  zero  valued  at  all  but 
one  point  on  the  time  axis.  Typical  values  for  the  impulse  are  1.2  lb-sec 
for  the  M16  Rifle  and  2.5  lb-sec  for  the  M79  and  M203  grenade  launchers. 

4.2  Factors  That  Affect  Weapon  Translation 

The  weapon  shoulder  support  is  represented  in  the  model  by  a linear 

dashpot  and  a piecewise  linear  combination  of  spring  elements,  which  results 

in  variable  shoulder  stiffness.  Thus  the  shoulder  spring  model  (illustrated 

in  Figure  18)  represents  the  effects  of  both  the  soft  compression  of  skin 

tissue  and  the  stretching  of  the  muscles  and  tendons  that  support  the  shoulder. 

The  shoulder  reaction  force,  F , is  an  explicit  function  of  coordinate  x, 

s 

but  since  the  man-weapon  system  is  dynamically  coupled,  the  shoulder  reaction 
force  also  is  influenced,  to  a lesser  extent,  by  both  the  pitching  coordinates 
0 and  <|>. 

Because  the  shooter  usually  holds  the  weapon  compressed  firmly  against 

his  shoulder  before  firing  and,  in  addition,  because  the  impulsive  breech 

force  compress  the  weapon  against  the  shoulder,  the  shoulder  support 

intuitively  is  expected  to  be  in  a state  of  strong  compression  (i.e.  x < x ) 

s 

for  most  of  the  firing  cycle.  This  hypothesis  was  substantiated  by  numerous 
computer  simulations  of  the  man-weapon  interaction  problem  and  also  by 
observations  of  experimental  load  cell  data  of  the  shoulder  reaction  force. 
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Figure  18  Shoulder  Spring-Dashpot  Model 


A flexible  shoulder  model,  consisting  of  a single  linear  shoulder 
spring  and  a dashpot,  therefore  provides  an  adequate  representation 
of  the  shoulder  impedance  to  weapon  recoil  loading. 

Experimental  results  also  indicate  that  the  shoulder  reaction  force, 
like  the  breech  force,  is  both  impulsive  and  nearly  periodic;  and,  in 
addition,  the  magnitude  of  the  peak  value  of  the  shoulder  force  tends  to 
decrease  for  successive  shots  in  a burst.  If  the  model  damping  coefficients 
are  near  critical,  then  the  shoulder  reaction  force  is  spread  over  a longer 
portion  of  the  firing  cycle,  and,  furthermore,  the  peak  values  of  the  shoulder 
reaction  force  become  sensitive  predominantly  to  the  damping  term  rather  than 
the  spring  stiffness  term.  In  order  to  duplicate  the  impulsive  and  the 
periodic  character  of  the  shoulder  force,  a damping  factor  between  10  and  30 
percent  of  critical  damping  is  required  (i.e.  c^  ■ 2 1 l/k^m^  where 
0.1  < < 0.3).  A spring  stiffness  of  250  lb/in  yields  pauk  values  for 

the  reaction  i'oKce  of  about  100  pounds.  Using  these  values  for  the  shoulder 
spring  and  dashpot  parameters , the  model  generated  shoulder  force  closely 
approximates  the  characteristics  of  the  experimentally  measured'  recoil  force. 

Results  of  the  computer  analysis  reveal  that  only  four  parameters  have 
any  significant  effect  on  the  duration  or  the  peuk  values  of  the  recoil 
force,  namely:  the  weapon  impulse,  I ; the  combined  mass  of  the  weapon  and 
the  arms,  m^.;  the  spring  stiffness,  kx;  and  the  damping  constant,  cx> 

Analog  computer  simulations  were  made  using  inertial  parameters  associated 
with  three  different  shooters  (small,  medium,  and  large)  and  these  results 
are  shown  in  Figures  19,  20,  and  21  respectively.  The  shoulder  dumping 
constant,  ranging  from  10  to  30  percent  of  critical  damping,  was  adjusted 
to  give  average  peak  recoil  forces  that  correspond  to  measured  results  from 
experimental  firings.  As  expected,  the  model  results  show  a decrease  in 
the  weapon  rotation  as  the  size  of  the  man  increases.  Only  slight  variations 
in  the  recoil  force  occur  among  the  three  cases. 

4 . 3 Factors  That  Affect  Weapon  Pitch  Rotation 

Certain  system  parameters  were  found  by  computer  analysis  to  have  a 
significant  effect  on  the  absolute  weapon  pitch  rotation,  while  other 
parameters  were  found  to  have  very  little  effect  on  weapon  rotation  at  all. 
Table  2 contains  a list  of  those  parameters  that  have  very  little  effect 
on  the  rifle  rotation.  The  torso  mass,  moment  of  inertia,  and  the  torsional 


TABLE  2 PARAMETERS  THAT  DO  NOT  AFFECT  WEAPON  ROTATION 


PARAMETER 

RIFLE 

ROTATION 

BEST 

APPROX. 

RIFLE 

ROTATION 

INCREASED 

VALUE 

RIFLE 

ROTATION 

kQ  (in-lb) 

2. 

3.84“ 

20. 

3.84° 

200. 

3.63" 

r C0 

0 cm.t 

0. 

3.84° 

0.5 

3.84“ 

1. 

3.83“ 

1cin<'lb“in_sec2') 

12. 

3.87“ 

26. 

3.84° 

50. 

3.83“ 

M(lb-eec2/ in) 

0.15 

3.77° 

0.3 

3.84“ 

0.45 

3.84“ 

L (in) 

10. 

4.07“ 

21. 

3.84* 

30. 

3.60“ 

rcm(ln> 

5, 

' 3.73 

13. 

3.84“ 

25. 

3.90* 

stiffness  and  damping  coefficients  for  the  8 coordinate  do  not  affect  the 
rifle  rotation  to  any  appreciable  amount  (a  few  hundreths  of  a degree  at 
most) . The  weapon  rotation  is  rather  insensitive  also  to  the  length  of  the 
torso  and  the  distance  from  the  torso  pivot  point  to  the  center  of  mass  of 
the  torso.  These  results  provide  some  insight  for  designing  a mount.  They 
suggest  that  a representation  of  the  torso  is  not  really  necessary  for  the 
design,  in  order  to  simulate  the  external  dynamical  motions  of  the  weapon. 

In  the  model,  the  motion  of  the  torso  mass  can  be  reduced  to  near  zero  by 
simply  assigning  a very  large  mass  and  moment  of  inertia.  For  example, 
by  increasing  the  mass  and  moment  of  inertia  by  a factor  of  one  thousand, 
the  rotation  of  the  weapon  increases  from  a value  of  3.84  degrees  to  a 
value  of  4.06  degrees  (see  Figure  22).  Although  the  torso  rotation  is 
reduced  from  7.0  degrees  in  the  nominal  case  to  zero  by  increasing  the 
mass  and  moment  of  inertia,  the  weapon  rotation  changes  only  by  0.2  degrees; 
consequently,  the  torso  mass  has  very  little  effect  on  the  total  rotation 
of  the  weapon. 

Figure  23  contains  a plot  of  the  weapon  pitch  versus  the  weapon  impulse 
after  0,3  seconds  have  elapsed.  Note  the  linear  relationship  between  these 
two  parameters  and  the  significant  effect  the  impulse  has  on  the  weapon 
rotation.  The  point  marked  "best  approximation"  on  this  and  on  other 
figures  corresponds  to  the  set  of  parameters  that  best  represents  the 
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dynamical  characteristics  of  the  165  pound  (medium  sized)  man.  Figure  24 
contains  a plot  of  the  weapon  rotation  at  0.3  seconds  and  at  0.15  seconds 
versus  the  firing  frequency.  Again  there  is  a linear  relationship  between 
the  weapon  rotation  and  the  firing  frequency.  Also  the  important  conclusion 
can  be  made  that,  for  a fixed  number  of  rounds  in.  a burst,  increasing  the 
firing  frequency  decreases  the  rifle  rotation  and  hence  lessens  the  weapon 
dispersion.  Figure  25  shows  the  effect  of  the  eccentricity  of  the  breech 
force,  6,  on  the  weapon  rotation.  This  parameter  represents  the  distance 
between  the  centerline  of  the  barrel  and  the  center  of  rotation  of  the  wea- 
pon at  the  shoulder.  The  weapon  pitch  rotation  varies  linearly  with  6. 

When  6 is  zero  no  torque  is  produced  about  the  shoulder  pivot  point;  never- 
theless, 1.5  degrees  of  rotation  occur  after  0.3  seconds,  which  is  due  mainly 
to  the  rearward  pitching  of  the  torso.  As  the  torso  pitches  rearward,  the 
rifle  pitch  angle,  <j>,  becomes  negative,  since  the  total  system  angular 
momentum  is  conserved  after  the  application  of  the  impulsive  breech  force. 

By  appropriately  positioning  the  stock  with  respect  to  the  barrel  centerline, 
the  rotation  of  the  weapon  can  be  virtually  reduced  to  zero. 

Parameters  that  have  a moderate  influence  on  the  rotation  of  the  wea- 
pon are  the  inertial,  stiffness  and  damping  parameters  associated  with  the 
man's  arms  and  the  weapon,  namely,  I , m , k,  , c , k , and  c . In  Figure  26 
the  rifle  rotation  is  plotted  as  a function  of  i/*^I  . A straight  line 
approximately  fits  the  points.  Similarly,  a plot  of  the  rifle  rotation 
as  a function  of "Vk^ , (see  Figure  27),  results  in  a linear  relationship. 

Thus  the  total  rifle  rotation  varies  linearly  with  both  l/^/l^  and*V . 

In  Table  3 the  rifle  rotation  at  0.3  seconds  is  1isted  foi  a range  of 

values  of  the  parameters  m , c,,  c , and  k . 

r’  x’  x 

Figures  28,  29,  30,  and  31  contain  the  output  of  a digital 
computer  simulation.  The  case  represented  is  that  of  a 165  pound  man 
firing  the  M16  rifle.  The  total  rotation  of  the  weapon,  shown  in  Figure  28 
is  the  sum  of  the  torso  rotation,  0,  and  the  weapon  rotation  relative  to 
the  torso,  (j>.  A characteristic  feature  of  the  simulation  is  the  jump  in 
acceleration  in  each  of  the  coordinates  after  each  shot,  due  to  the  impulsive 
breech  force. 
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TABLE  3 SENSITIVITIES  OF  PARAMETERS  THAT  AFFECT  WEAPON  ROTATION 


PARAMETER 

INCREASED 

VALUE 

RIFLE 

ROTATION 

BEST 

APPROX. 

RIFLE 

ROTATION 

DECREASED 

VALUE 

RIFLE 

ROTATION 

in^  (lb-scc^/ln) 

0.0!) 

4.48° 

O.0S9 

3.84“ 

0.2 

2..' 

crit 

0. 

3.69° 

.5 

3.84“ 

1. 

4.02 

ks (lb/in) 

5. 

5.07“ 

23. 

3.84“ 

250. 

2.99 

co^ca  crit 

0. 

5.29“ 

1. 

3.84” 

1.5 

3.55 

4.4  Recommended  Design  Configurations 

Based  on  experimental  findings  and  model  predictions,  two  possible 
design  configurations  for  a universal  small  arms  mount  are  recommended 
and  are  illustrated  in  Figures  32  and  33.  By  appropriate  choices 
for  the  masses,  lengths,  viscous  damping  constants,  and  spring  stiffnesses, 
the  mount  can  be  designed  such  that  its  dynamical  characteristics  are 
analogous  to  the  characteristics  of  the  simulation  model.  The  simplified 
version  of  the  design  in  Figure  33  is  based  on  the  fact  that  the  rotation 
of  the  torso  mass  does  not  significantly  affect  the  weapon  motion  or  the 
shoulder  reaction  force. 
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5 . 0  CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  Conclusions 

a.  Based  on  experimental  results  and  model  predictions,  the  force 
reaction  at  the  shoulder  can  be  approximated  by  a single  linear  spring 

and  dashpot.  By  the  proper  adjustments  of  the  spring  and  damper  constants, 
both  the  impulsive  and  periodic  behavior  of  the  reaction  force  can  be 
approximated  as  well  as  the  magnitudes  of  the  peak  values. 

b.  The  parameters  that  affect  the  peak  forces  at  the  shoulder  are 

the  shoulder  spring  stiffness,  the  shoulder  damping  coefficient,  the  impulse 
of  the  weapon,  and  the  effective  mass  of  the  weapon,  arms,  and  hands. 

Other  model  parameters  have  only  a minor  influence  on  the  peak  values  of 
the  shoulder  reaction. 

c.  The  weapon  rotation  after  a specified  time  is  not  sensitive  to 
the  torso  inertial  or  geometric  parameters.  By  effectively  removing  this 
link  in  the  model,  the  weapon  rotation  is  not  significantly  changed.  The 
parameters  that  affect  the  weapon  rotation  most  are  the  weapon  impulse, 
the  eccentricity  of  the  barrel  centerline  from  the  shoulder  pivot  point 
and  the  firing  frequency.  The  weapon  rotation  was  found  to  vary  linearly 
with  each  of  these  parameters.  Parameters  that  affect  the  weapon  rotation 
but  are  of  secondary  significance  are  the  effective  inertial  parameters  for 
the  weapon  and  arms,  the  torsional  stiffness  and  damping  of  the  weapon 
pitch,  and  the  shoulder  stiffness  and  damping  constants. 

d.  The  recommended  configuration  for  a mount,  based  on  experimental 
and  analytical  results,  is  presented  in  Figure  33.  It  allows  for  two 
degrees  of  freedom  (pitch  and  translational)  and  contains  a spring  and 
damper  to  simulate  the  shoulder  reaction  force  and  an  effective  torsional 
spring  and  damper  to  simulate  the  man's  torsional  stiffness  and  damping 
characteristics . 

5.2  Recommendations  For  Future  Work 

The  recommendations  for  future  work  in  developing  a universal  small 
arms  mount  design  are  presented  below: 

a.  Additional  tests  should  be  performed  to  support  the  complete  design 
of  a universal  small  arms  mount  fixture  or  the  modification  of  an  existing 
design.  The  tests  would  be  performed  to  refine  the  values  obtained  for 
the  design  parameters. 
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b.  If  further  tests  are  conducted,  a better  method  for  measuring  the 
shoulder  reaction  force  would  be  to  cut  off  the  end  section  of  weapon  stock 
and  insert  one  or  two  force  washers  and  an  accelerometer  directly  to  the 
stock.  This  method  would  replace  the  need  for  the  rather  heavy  and 
cumbersome  shoulder  force  measuring  device  shown  in  Figure  15  and  would 
therefore  result  in  less  errors  in  measurement. 

c.  Having  additional  experimental  data,  iefine  the  mount  parameters, 
using  the  biodynamics  model  to  best  approximate  the  actual  shooter's 
response.  Then  incorporate  the  refined  parameter  values  in  the  mount  design. 

d.  Investigate  existing  mount  designs  to  determine  if  modifications 
can  be  made  on  them  to  satisfy  the  design  requirements  of  the  recommended 
concept.  If  an  existing  mount  is  found  that  can  be  modified  ro  meet  these 
requirements,  then  the  authors  recommend  that  it  be  adapted  to  the  design 
configuration  in  Figure  33  and  tested. 

e.  Otherwise,  if  no  existing  design  can  be  so  modified,  then  a 
universal  small  arms  mount  should  be  designed  from  scratch  based  on  the 
concept  shown  in  Figure  33  and  tested. 
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APPENDIX  A 

A.O  MODEL  REPRESENTATION  OF  THE  HUMAN  BODY  SEGMENTS 
A. 1 Head  and  Torso  Representation 

The  formulas  used  to  calculate  the  human  body  segment  parameters  and 
the  methods  for  combining  the  body  segments  to  represent  the  torso  and  arms 
are  presented  below.  The  head  and  trunk  model  representation  is  shown  in 
Figure  A-l.  An  ellipsoid  of  revolution  represents  the  head,  while  ellip- 
tical cylinders  represent  the  trunk.  Each  of  the  distances  illustrated  in 
Figure  A-l  are  defined  by  certain  anthropometric  dimensions.  In  particular, 

a.  d = \ (stature  - cervical  height) 

bj^  = (head  circumference)  /2tt, 

at  » jr  (chest  breadth  + waist  breadth  + hip  breadth) 
b = (chest  depth  + waist  depth  + hip  depth) 

£ut  = cervical  height  - substernal  height, 

and 

m substernal  height  - trochanteric  height. 

The  mass  and  moments  of  inertia  for  the  head  are  given  by  the  standard 
formulas  for  an  ellipsoid  of  revolution, 

= 3 ?r  ahd  bhd  phd’ 

^aahd  ^ ^hd^  * 

and 

Ibb,  , “ IccLJ  “ “hd  (ahd  + bLP  /5 
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b 
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Figure  A-l  Head  and  Torso  Rigid  Body 
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and  the  mass  density 


= 1.11  x 103 


(kg/m3) 

was  obtained  from  reference6.  For  the  upper  torso  the  mass  and  moments 


of  inertia  are  given  by 


and 


m = it  a.  b Z p 
ut  t t ut  ut 


l-  = mut  (at  + bt}  74 


aa 


ut 


hb  ■ "ut  ° "1  + *it>  nz 
ut 


I = m „ (3  b2  + £2  ) / 12 
cc  ut  v t ut 

ut 


Again  the  mass  density 

Put  * 9.2  x 102  ^kg/m3) 

was  obtained  from  reference6,  page  195.  Likewise  for  the  lower  torso 
m£t  = * at  bt  llt  p£t 


Xaa£t  = mlt  (a?  + b^  /4 
Ibb£t  “ mlt  ° at  + 712 


= “et  <3  bt  + 712 


-it 


and 


p£t 


- 1.01  x 10; 


(kg/m3) 


In  the  model  tue  head,  upper  torso  and  lower  torso  are  combined  into  a 
single  rigid  body  with  mass  M, 

M ’ "hd  + "ut  + mlt  ■ 


®Dempstcr,  W.T,  , "Space  Requirements  of  the  Seated  Operator,"  Wright 
Air  Development  Center,  TK-55-159,  Wriglit-Pattersori  Air  Force  Base, 
Ohio,  AD  78792,  1955 


Furthermore,  the  center  of 
rcm  • lmhd  Uit  * At 


mass  location,  rcm,  is  given  by 
+ ahd>  + mut  (££t  + °‘5  £ut) 


+ °*5  m£t  ££t]  /M 

and  the  moment  of  Inertia,  I , about  the  center  of  mass  of  the  head 

cm 

and  trunk  rigid  body  is  given  by, 


cm 


= I 


bb 


+ I 


hd 


bb 


+ I 


ut 


bb 


'It 


+ m,  , (r  - £ . - £«.  - a,  + m . (r  - £p. 
nd  cm  ut  £t  hd'  ut  v cm  £ t 

+ m n (r  - 0.5  £ )2 

£t  cm  ut 


0.5  lut> 


A. 2 Right  Arm,  Left  Arm  and  Weapon  Configuration 

The  geometrical  representations  for  the  upper  arm  segments,  forearm 
segments  and  hands  are  shown  in  Figure  A-2.  The  hand  is  represented  by 
a sphere  with  density 

“ 1.148  x 103  (kg/m3) 

obtained  from  reference9  and  mass, 

m.  = 0.01  W + 0.1588  (kg) 

where  W = total  body  weight,  obtained  from  reference7. 

The  radius  and  moment  of  inertia  formulas  for  the  spherical  hand  are 
(3  n^)1^ 
k 3/4  u 

and 

jn  = J22  = J33  - f ^ 

h h h 

A frustum  of  a right  circular  cone  is  used  to  represent  both  the 
upper  arm  segment  and  the  forearm  segment.  Lenghts  £u  and  £^  are  defined 
by  anthropometric  dimensions, 

£u  = acromial  height  - radiale  height 


9Drillis,  R.  and  Contini,  R.  , "Body  Segment  Parameters,"  Office  of  Voca- 
tional Rehabilitation,  Department  of  Health,  Education  & Welfare,  N.Y. 
University,  School  of  Engineering  & Scienre,  N.Y.,  Report  1166-03,  1966 

7Barter , J.T. , "Estimation  of  the  Hass  of  Body  Segments,"  Wright  Air 
Development  Center,  Wright-Patterson  Air  Force  Base,  Ohio,  TR-57-260, 

AD  118222,  1957 
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Figure  A-2  Geometrical  Arm  and  Hand  Representations 
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and 


£_£  = radiale  height  - stylion  height 
coefficients  p^  and  were  obtained  from  reference^; 


Pu  =*  0.449 


P L = 0.423 


and  the  masses  and  m ^ of  the  upper  arm  and  forearm,  respectively  were 
obtained  from  reference7: 


mu  = 0.04  W - 0.6577 


(kg) 


and 


m£  = 0.02  W - 0.1134 


(kg) 


where  W " total  body  weight  (kg) . 

The  moments  of  inertia  of  the  upper  arm  are  then  given  by  the  standard 
formulas  for  a frustrum  of  a cone, 


aa 


2 A m2 
u u 

p T~ 
a u 


and 


bb 


J = m [B  lz  + A m /(p  l ] 
cc^  u u u u u a u 


where 


u 


+ P3  + p4 
a a a 

20  tr  ic2 
a 

3 [Ka  + 3ya  Ka  + (5  + Ka>  ^ 

80  k2 
a 


and 


4 p -1 
u 


1-2  Pu  + /12  pu  (1  - pu)  -2 


1 + p + p2 

a a 
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Moments  of  inertia  formulas  for  the  forearm  are  identifical  to  those  of 
the  upper  arm  except  that  the  coefficients  Z , Pu  and  m^  are  replaced 
and  m£  respectively. 

Figure  A-3  contains  an  illustration  of  the  two-dimensional  geometrical 
configuration  of  the  weapon  and  both  arms.  Point  S locates  the  shoulder 
pivot  point.  For  convenience  the  assumption  is  made  that  both  arms  are 
connected  at  point  S.  The  right  hand  grasps  the  weapon  at  distance  ci 
and  the  left  hand  grasps  the  weapon  at  distance  C2  and  distance  locates 
the  center  mass  of  the  weapon.  Table  A-l  contains  the  weapon  mass,  moment 
of  inertia  and  center  of  gravity  location  of  the  M16  rifle  and  the  M203 
grenade  launcher,  as  measured  by  Mr.  T.  Hutchings. 


TABLE  A-l  WEAPON  PARAMETERS 


PARAMETER 

M-16 

M2  03 

m (lb) 

w 

8.6 

9.4 

I (lb-in-sec2) 
w 

1.6 

3.6 

X (in) 

w 

17. 

i 

19. 

Referring  to  Figure  A-3,  the  total  mass  of  the  arm-weapon  system  and  the 
center  of  mass  location  (xq,  yQ)  are  given  by 

m = 2 m + 2mn  + 2m,+m 

r u Z ti  w 

x = [m  p Z (cos  a + cos  3) 
o u u u ' 

Lg)  (cos  a + cos  3) 

+ “h  (ci  + C2)  + mw  x]  /mr 

p£  (sin  a + sin  ^ /mr 


[m  p Z + m„  (Z  + 
u ru  u Z u 


where 


cos  3 = + c^  - 

2 Z ci 

u 1 


" £U  + C2  ~ 

1 lu  c2 


ml  + PZ 


with  Z^t  p£ 


cos  a 


Finally  the  combined  moment  of  inertia  of  the  arm-weapon  system  about  the 
center  of  mass  becomes 


Jcm  ' 2 Jn  + 2 Jbb  + 2 Jbb.  + J 

h u JL 


w 


+ [<xo  - Cl)2  + (xQ  - c2)2  + y2] 

+ mu  [(xo  “ pu  £u  C0S  a)2  + (xo  ' pu  £u  cos  3)2 

+(yo  “ pu  £u  Sin  a)2  + (yo  " pu  £u  Sin  3)2  ] 

+ m^  {[xQ  - (£u  + p£  cos  a]2  + [xq  - Uu  + p£  lg)  cos  0]2 


+ [yQ  - Uu  - P£  lg)  sin  a]2  + [yQ  - - p ^ £^>)  sin  g]2j 

+ mw  t(xo  " \,)2  + yo] 


A. 3 Anthropometric  Data  and  Calculated  Results 

A list  is  provided,  in  Table  A-2,  of  some  typical  weights  and  anthro- 
pometric dimensions  for  three  men  of  various  sizes  (small,  medium,  large). 
This  table  contains  a partial  list  of  data  that  were  measured  by  Santschi 
reference8.  Table  A-3  contains  the  calculated  dimensions  and  inertial 
parameters  based  on  the  data  in  Table  A-l  and  A-2.  The  arm-weapon  inertial 

parameter,  I , is  somewhat  larger  than  actual  because  the  five  pound 
cm 

shoulder  force  measuring  device  was  included  in  the  inertia  and  center  of 
mass  calculations.  The  results  in  Table  A-3  show  the  orders  of  magnitude 
of  the  model  parameters  along  with  the  expected  variations  for  different 
size  men. 


^Santschi,  W.R. , DuBois,  J.,  and  Omoto,  C. , "Moments  of  Inertia  and 
Centers  of  Gravity  of  the  Living  Human  Body,  Aerospace  Medical 
Research  Laboratories,  Wright-l’atterson  Air  Force  Base,  Ohio, 
AMRL-TDR-63-66,  AD  410451,  1963 
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TABLE  A- 2 (SANTS CHI) 

WEIGHT  AND  ANTHROPOMETRIC  DATA  FOR  THREE  TEST  SUBJECTS 


Weight 
Stature 
Head  Circ. 

Chest  Breadth 
Waisr.  Breadth 
Hip  Breadth 
Chest  Depth 
Waist  Depth 
Hip  Depth 
Cervical  ht. 
Substernal  ht. 
Trochanteric  ht. 
Forearm  Length 
Upper  Arm  Length 


SUBJECT  1 


161.6  (lb) 
66.9  (in) 
22.2 
13.3 
11.6 
13.8 
9.3 


SUBJECT  2 


172.2  (lb) 

73.0  (in) 

23.5 

12.7 

12.0 

14.8 
10.3 

9.3 

10.6 

62.7 

50.8 

39.2 
11.5 

13.7  NF 


SUBJECT  3 


210.4  (lb) 

72.3  (in) 

23.3 

15.6 

13.2 

16.3 

10.4 

10.9 

10.5 

62.4 

50.7 

38.5 
11.1 
14.0  v 


TABLE  A- 3 CALCULATED  BIOMECHANICAL  PARAMETERS 
FOR  THREE  TEST  SUBJECTS 


I (lb-in-sec2) 
cm 


1^,  (lb- in-sec2) 
L (in) 
ron  (ln) 

*„  <ln> 


SUBJECT  //I 

SUBJECT  #2 

] 

SUBJECT  if 3 

87. 

97. 

116. 

18. 

21. 

26. 

32. 

33. 

38. 

9.7 

11. 

13. 

20. 

21. 

21. 

13. 

14. 

13. 

9.6 

9.5 

9.3 

ooo  noon  non  non  non  non  o o O’  oo 


APPENDIX  B 


B.O  DIGITAL  COMPUTER  PROGRAM  LISTING 


REAL  M.MR.ICM.IR.L.  KXW ,KXS  .KT.KP.MU.NU.K.IF 

COMMON  /$A/  M.MR.ICM.IR.L.O.XS.RCM.KXW.KXS.KT.KP.CDX. 

* COT.CDP.G.K.XO.TO.PO.STO.STPO.CTPO.TIM 
COMMON  /$PL0T/  NP . TIMEC 502 ) ,R( 502 ) , Z( S02 ) 

0IMENSI0N  PRMT(S)  ,RHS(6) . YC  6 3 . AUX( 8 .6 ) .ERW( 6 ) 

EXTERNAL  FCT.OUTP 
CALL  PLOTS  ( IBUP . I0UM.14) 

ALL  UNITS  ARE  IN  THE  POUND.  INCH.  SECOND  SYSTEM 

G = 386.4 
NDIM  - S 
NRUN  = 0 

READ  INERTIA  OATH 

10  READ  ( S . 1000 ,END=999 3 M.MR.ICM.IR 
NP  = 0 

NRUN  = NRUN  + 1 

READ  DISTANCES 

READ  (S.1000)  L.D.XS.RCM 

READ  SPRING  CONSTANTS 

READ  (5.1000)  KXW.KXS.KT.KP 

READ  DAMPING  RATIOS 

READ  (5.10001  CRX.CRT.CRP 

CALCULATE  DAMPING  COEFFICIENTS 

CDX  = CRX*2.wSQRTCKXS*MR) 

COT  = CRT*Z.*SQRT(KTwICM) 

CDP  = CRPw2.*SQRT(KP*IR) 

READ  INTEGRATION  PARAMETERS:  START  TIME.  STOP  TIME. 

TIME  INTERVAL.  UPPER  ERROR  BOUND.  AND  PRINT  PARAMETER 

READ  (5.1000)  ( PRMT (I).I— 1.4) 

READ  INITIAL  CONDITIONS:  XD0T(l-3),  X£ 1-3 ) 
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nonnnn  non  non  non 


C 

C 

C 


READ  (5.1 000  ) Y 

READ  ERROR  WEIGHTS  FOR  RKOS 
READ  (5,1000)  ERW 

READ  SHOT  DURATION.  AND  THE  TIME  BETWEEN  SHOTS  AND  WEAPON  IMPULSE 

READ  (5.1000)  TOUR. TINT. IF 

PRINT  INPUT  DATA 

WRITE  ( 6,2000 ) NRUN 
WRITE (6, 3000)  M.MR.ICM.IR 
WR I TE( 6 , 4000  ) 

WRI TE( 6,3000  ) L.D.XS.RCM 
WR I TE( 6,5000 ) 

WRITE( 6 ,3000  ) KXW ,KXS  ,KT  ,KP 
WR I TE( 6,6000  ) 

WR I TE ( 6 . 3000  ) CRX.CRT.CRP 
WRITE( 6 ,7000 ) 

WRITE( 6 .3000  ) ( PRMT( I ) . 1 = 1 .4 ) 

WR I TE ( 6 ,8000 ) 

WRI TE ( 6.3000 ) ( YC  I ) , 1 = 1 ,3) 

WR I TE ( 6 ,9000  ) 

WR I TE ( 6 , 3000 ) (Y(I),I=4.6) 

WR I TE( 6,5500 ) 

WRITE( 6 ,3000 ) ( ERW (I) . I = 1 . 3 ) 

WR I TE( 6 , 6500 ) 

WRI TE( 6 .3000  ) ( ERW( I ) . 1=4 .6 ) 

WRI TE ( B ,7500 ) 

WRITE ( 6.3000 ) TDUR , TINT . I F 
WRITE  HEPOINO  FOR  OUTP 


1 

i 


WR I TE ( 6 , 1 600  ) 

INITIALIZE  PARAMETERS 

TEND  = PRMT ( 2 ) 

K = ICM  + MwRCMwmZ 

CALCULATE  INITIAL  CONDITIONS  AFTER  THE  NEXT  SHOT 


ilii  iriilrl 


c 

c 

c 


20 


25 


SO 


60 


C 

C 

C 


xo  = 

Y(  4) 

TO  = 

YC  5 ) 

PO  = 

Y ( 6 ] 

STO 

= S INC  TO ) 

CTO 

= COSC  T03 

SPO 

= SlNCPO) 

CPO 

= cost  PO  3 

STPO 

= SINCTO 

CTPO 

= COSCTO 

MU  = 

MRwLmCPO 

DL  = 

LmCPO  + 

CALL 

OUTPCT.Y 

PRMT ( 1 J = TIM 

PRMTC  2 ) = TIM 

PO] 

PO) 


+ TDUR 
+ TINT 

TEND]  00  TO  50 
TEND]  PRMTC2)  = TEND 


IFCPRMTC 1 ) .OE 
IF(  PRMT ( 2 ) .OT 
DO  25  1=1.3 
RHSU)  = ERWCI) 

ETfl  = Y ( 4 1 + LkSPO 
NU  = IR  + MR*Y( 4)*w2 

Bll  = K»NU  MRwIRwLi«*2  + ( MU*Y(  4 ) )»»2 

812  = NUkMU 

813  a -(IR  + MRwYC  4)*ETR)wMU 
B22  = MRwNl) 

823  a -MR» ( I R + MRwY ( 4 )*ETA ) 

B33  = MRwCK  + IR  + MR»ETA>m»2  ] 

DELTA  = KwMRkNLJ  + IRa(MR»LiiSP0]»n»2 
Y ( 1 ) = Y( 1 ] + IFw ( -8 1 1 + 812wDL  + B 13*D  J/DELTA 

Y(  2 ) = Y(  2 ] * IFw(-B12  + B22*DL  + B23*D 3/DELTA 

Y ( 3 ) = Y( 3 ) + IFw( -B13  + B23««DL  + B33*D)/DELTA 

INTEGRATE  THE  MAN/WEAPON  EQUATIONS 

CALL  RKGSC  PRMT . Y ,RHS .NDIM . I HLF ,FCT .OUTP ,AUX ) 
WRITEC  6 .2500]  I HLF 
00  TO  20 
WRI TEC  6 .3500 ] 

00  60  1=1 ,NP 

WRI TEC 6. 4500]  I .TIMEC 1 3 ,R( I ) ,Z( I ) 

CONTINUE 

PLOT  RESULTS 

CALL  CALPLT ( NRUN  3 
00  TO  10 
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1 1 1'ihiTur  mu  ii 


999  CALL  PLOT( 12.0.0.0.999) 

1000  FORMRTC 8F10 .0 ) 

1500  FORMAT!  * 1 \T4 . 'TIME ' .T17 , 'XOOT  * .T37 . ‘TOOT  ’ .TS7 , 'PDOT  * .T77.  'X  ' , 

* T97 , 'T' .T117.  ’P') 

2000  FORMAT!  *1  THREE  DEGREE  OF  FREEDOM  MAN/WEAPON  MODEL  INPUT  * .5X . 'RUN  ' 

* .13/ 

* 'OINERTIA  DATA: */T3» 'M*.T2B. ’MR'.T53, * ICM  * , T78  » " I R ’ J 
2600  FORMAT!  'OIHLF  =',I3//) 

3000  FORMAT! IX . 4! E IS  .8 . 1 OX  ) ) 

3600  FORMAT!  *1  NP '. T10 . 'TIME T30 , 'RECOIL '. TSO . 'PI TCH * ) 

4000  FORMAT!  '00 1 STANCES! */T3 . ‘L ' .T28 . 'D  * , T53 . *XS * , T78 . ’RCM ’ ) 

4500  FORMAT!  ’ ‘ . 13 . T 10 , F8 .5 . T30 ,E IS . 8 . T50 ,E 1 5 .8  ) 

5000  FORMAT!  'OSPRINO  CONSTANTS! */T3, 'KXW  '.T28. 'KXS  ' . T53 . 'KT * . T78  . *KP  ' ) 
5600  FORMAT!  'OVELOCITY  ERROR  WEIGHTS : */T3 . 'XOOT  '. T28 . 'TOOT *.  T53 , 'POOT  ' ) 
6000  FORMAT!  'OCRITICAL  DAMPING  RATIOS : */T3 . ‘CRX  T28 . 'CRT  \ T53 . 'CRP * ) 
6500  FORMAT!  'OPOSITION  ERROR  WEIGHTS: '/T3. 'X'.T28.  ’T'.T53, *P  * ) 

7000  FORMAT!  'ORKGS  PARAMETERS : ’/T3 , 'START  TIME '. T28 . 'STOP  TIME’. 

* T53 . 'TIME  INCR  ' ,T78. 'ERROR  BOUND') 

7600  FORMAT!  'OTIME  PARAMETERS! '/T3. 'TOUR  ' , T28 . ' TINT  ' . T53 . ' IF  ’ ) 

6000  FORMAT!  'OINITIAL  VELOCITIES: */T3. 'XDOT ' .T28.  'TOOT  ' , TS3 . 'PDOT  ' ) 

9000  FORMAT!  'OINITIAL  POSITIONS:  '/T3,  *X'.T<c8.  ' T ' *T53 , 'P  ' ) 

STOP 

END 
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SUBROUTINE  CALPLTCNRUN) 

COMMON  /$PL0T/  NP , TlMEf 502  ) ,RC S02 ) ,Z( S02 ) 

DIMENSION  TC4).DAT(7) 

Til)  = TIME! 1 ) 

TC2)  = TIME(NP) 

FRN  = NRUN 
C 

C PREPARE  FOR  THE  FIRST  PLOT 

C 

CALL  NEWPEN(l) 

CALL  PLOTCO. 0,-15. 0.-3) 

CALL  DASHPCO.0.16.0.0.3) 

CALL  PLOTC 1 .0. -IS. 0.-3) 

CALL  PLOTCO. 0.3. 0,-3) 

CALL  RECT(0.0,0.0,7.5,9.0,a.0.2) 

CALL  PLOTC 1 .5,0.76.-31 
CALL  SCALEC  T .7 .0 ,2 . 1 ) 

CALL  SCALEC  R ,5 .26 ,NP , 1 ) 

CALL  SCALEC  Z ,5 .25 ,NP . 1 ) 

TIMEC  NP  + 1 ) = TO) 

TIMEC  NP+2 ) = TC  4- ) 

C 

DRAW  AXES  AND  GRAPHS 
C 

CALL  AXISCO. 0.0.0. 'TIME  ( SECONDS  .1  • . -1 4 .7 .0.0 .0  . TC  S ) . TC  A ) .* 

CALL  AXISCO. 0,0.0, 'PITCH  ANGLE  C DEGREES ) ’ .21  ,5 .26 .90 .0 . Z( NP  + n 
* Z(NP*2)) 

ZEROP  = ABSC ZC  NP+1  )/Z( NP+2  3 ) 

CALL  PLOTCO. 0. ZEROP. 3) 

CALL  PLOT. '7.0, ZEROP. 2) 

CALL  LINECTIME.Z.NP.l  ,0.1 ) 

CALL  NEWPENCE) 

CALL  AXISC-0. 76, 0.0. 'MOUNT  FORCE  C POUNDS  ) ’ ,20 . 5 .26 .90 .0 , 

•*  RCNP+l  ) ,R( NPv2 ) ) 

ZEROP  = ABSCRC  NP  + U/RC  NPf2  ) ) 

CALL  PL0TC-0.76.ZEROP.3) 

CALL  PLOTC  7 .00 .ZEROP ,2 ) 

CALL  LI  NEC  TIME .R  ,NP . 1 ,0.1) 

C 

C WRITE  TITLE  HEADINGS 

C 

CALL  NCWPENC 1 ) 

CALL  SYMB0LC1 .5,6.26,0.21 . ‘MAN/WERPON  ANALYSIS  * .0.0. 19) 

CALL  SYMBOLC 1 .5,6.0.0.14. 'MOUNT  FORCE  AND  TOTAL  PITCH  ANOLE ' . 

» 0.0.33) 


• : -5 
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CALL  OBDATE!  DAT ) 

CALL  SYM80LC 1 .5, 5. 75, 0. 14.  DAT. Q. 0.28) 
CALL  SYMBOL! 1 .5.5.50.0.14. ‘RUN  ‘,0.0,4) 
CALL  NUMBER! 999. .999. ,0.14. FRN. 0.0.0) 

PREPARE  FOR  THE  NEXT  PLOT 

CALL  PLOT (8. 5, -15. 0.-3) 

CALL  DASHP (0.0, 16-0. 0.3) 

RETURN 

ENC 


SUBROUTINE  FCTIT.Y.RHS) 

REAL  M.MR.ICn.IR.L.KXW.KXS.KT.KP.nU,NU,K 

COMMON  /*R/  M .MR . ICM . IR .L .0 ,XS .RCM.KXW ,KXS *KT ,KP ,COX . 

* CDT.CDP.C.K.XO.TO.PG.STO.STPO.CTPO.TIM 
DIMENSION  YI6J .XI3) .XDI3) .RHSCS) 

DO  10  1=1.3 
J=I+3 

XU)  = Y(J) 

10  XDII)  = Y ( I ) 

STP  = SINI X( 2)  + X£31) 

C-TP  = COS(  XI 2 ) ■*•  Xt  3 ) 3 
ST  = SINI XI 2 ) ) 

CT  = COSI XC2  ) ) 

SP  = SINI XI 3 X ' 

CP  = CCS! XI 3 ) ) 

MU  = MR*L*CP 

ETfl  = XI  1 ) •*•  L*SP 

NU  = IR  + MRwXC 1 )**Z 

611  = K*NU  i-  MR*IR*L**2  + I MU*X(  1 ) )w*2 
B12  = NUmMU 

B13  = -CIR  MR»X(  1 )*ETfl )*MU 

822  = MRwNU 

823  = -MRwIIR  + MR»X( 1 )«ETfl) 

833  = MR*(K  + IR  + MR»ETH«<*Z ) 

DELTA  = K*MR*NU  + IR*(  MR*L*SP  ) 

EX  = MR»X(  1 )*( XOI 2 ) * XDODwxZ  + MR*L*SP»XD(Z)**2  +MR*D*(STPO 
p -STP) 

FKX  = KXWalXin  - XO)  + COX  wXD(  1 ) 

IFIXt 1) .LT.XS)  FKX  = KXS*( X(1)-XS )+KXW»( XS-XO )+CDX  *XD( 1 ) 

FX  = FX  - FKX 

FT  = -2.*MR*ETfl*XD(  1 J*(  XOI 2 )+XD(  3 ) )-X(l  )*MU*XDC  3 )*(  2 .*XD(  2 )<-XD(  3 ) ) 
« -CDT*XD(2)-KT*(X(2)-T0)  *(  M*RCM+MR*L  )*0>»(  ST-STO  ) 

* -MR *0*1 XI 1 )*CTP-XO“CTPO) 

FP  = -2.*MR*X( 1 )*XDt l)*(XD(Z)+XD(3))vMU*XD' 1 )*XD( 2 )**Z-CDP*XD( 3 ) 

* -KP*( XI 3 )-P0 )-MR»0*X(  1 )*CTP  f MR*0*XO*CTPO 

RHSI 1 ) = ( FXwB 1 1 + FT*B 12  + FP*B  1 3 )/DELTfl 

RHSC2)  = (FX*B12  + FT*B22  + FP*B23  J/DELTR 

RHSI 3 ) = I FX*B13  + FT*B23  + FP*B33 J/DEITA 

RHSI 4)  = X D C 1 ) 

RHSI 5 ) = XDC  2 ) 

RHSI 6 ) = XOI 3 ) 

RETURN 

END 


SUBROUTINE  GUTP( T. Y ,RHS . IHLF .NDIM.PRMT I 

REAL  M.MR.ICM.IR.L.KXW.KXS.KT.KP.MU.NU.K 

COMMON  /$A/  M.MR.ICM.IR.L.O.XS.RCM.KXW.KXS.KT.KP.COX. 

« CDT.COP.G.K.XO.TO,PO,STO,STPO.CTPO,TIM 
COMMON  /$PLOT/  NP . TIME! 502 ) .RC 502 ) .1 ( 502 I 
DIMENSION  Y f 6 ) .RHSC  6 I ,PRMT( 5 I 
TIM  = T 

CPI  = 1S0./3.14159Z7 
NP  = NP  1 

IFCNP.GT.50DI  GO  TO  50 

RC  NP I = KXW* ( Y( 4 ) - XOJ  + CDX  »Y(1) 

IFCYC4)  .LT.XS)  RCNPI  = KXS»(YC4)  - X8)  + KXWwCXS  - XO ) + CDX  *Y{ 1) 
ZCNP)  = CYC5I  + Y(6))*CPI 
TIMECNP)  = T 
50  CONTINUE 

C ADDITION  FOR  OTHER  PLOTS  »»»»«»»»»» 

HR  I TEC  3 IT  . Y ,RHS 
WRITEC6 . 1000 1 T.Y 

1000  FORMAT ( 1X.F8.5.T10.6C5X.E15.8)) 

RETURN 

END 
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SUBROUTINE  RKGS 

PURPOSE 

TO  SOLVE  R SYSTEM  OF  FIRST  ORDER  ORDINARY  DIFFERENTIRL 

EQUATIONS  WITH  GIVEN  INITIAL  VALUES. 

USAGE 

CALL  RKGS  ( PRMT  ,Y .DERY ,NDIM . IHLF ,FCT . OUTP .AUX ) 

PARAMETERS  FCT  AND  OUTP  REQUIRE  AN  EXTERNAL  STATEMENT. 

DESCRIPTION  OF  PARAMETERS 

PRMT  - AN  INPUT  AND  OUTPUT  VECTOR  WITH  DIMENSION  GREATER 
OR  EQUAL  TO  6.  WHICH  SPECIFIES  THE  PARAMETERS  OF 
THE  INTERVAL  RND  OF  ACCURACY  AND  WHICH  SERVES  FOR 
COMMUNICATION  BETWEEN  OUTPUT  SUBROUTINE  (FURNISHED 
BY  THE  USER)  AND  SUBROUTINE  RKGS.  EXCEPT  PRMTCS) 
THE  COMPONENTS  ARE  NOT  DESTROYED  BY  SUBROUTINE 
RKGS  AND  THEY  ARE 

PRMTU)-  LOWER  BOUND  OF  THE  INTERVAL  (INPUT). 

PRMT( 2 )-  UPPER  BOUND  OF  THE  INTERVAL  (INPUT), 

PRMT( 3 )-  INITIAL  INCREMENT  OF  THE  INDEPENDENT  VARIABLE 
(INPUT). 

PRMT( 4 )-  UPPER  ERROR  BOUND  (INPUT).  IF  ABSOLUTE  ERROR  IS 
GREATER  THAN  PRMTU).  INCREMENT  GETS  HALVED. 

IF  INCREMENT  IS  LESS  THAN  PRMT(3)  AND  ABSOLUTE 
ERROR  LESS  THAN  PRMT(4)/SO.  INCREMENT  GETS  DOUBLED 
THE  USER  MAY  CHANGE  PRMTU)  BY  MEANS  OF  HIS 
OUTPUT  SUBROUTINE. 

PRMT(S)-  NO  INPUT  PARAMETER.  SUBROUTINE  RKGS  INITIALIZES 
PRMT(5)=0.  IF  THE  USER  WANTS  TO  TERMINATE 
SUBROUTINE  RKGS  AT  ANY  OUTPUT  POINT.  HE  HAS  TO 
CHANOE  PRMT(S)  TO  NON-ZERO  BY  MEANS  OF  SUBROUTINE 
OUTP.  FURTHER  COMPONENTS  OF  VECTOR  PRMT  ARE 
FEASIBLE  IF  ITS  DIMENSION  IS  DEFINED  GREATER 
THAN  6.  HOWEVER  SUBROUTINE  RKGS  DOES  NOT  REQUIRE 
AND  CHANOE  THEM.  NEVERTHELESS  THEY  MAY  BE  USEFUL 
FOR  HANDING  RESULT  VALUES  TO  THE  MAIN  PROORRM 
(CALLING  RKOS ) WHICH  ARE  OBTAINED  BY  SPECIAL 
MANIPULATIONS  WITH  OUTPUT  DATA  IN  SUBROUTINE  OUTP. 

Y - INPUT  VECTOR  OF  INITIAL  VALUES.  (DESTROYED) 

LATERON  Y IS  THE  RESULTING  VECTOR  OF  DEPENDENT 
VARIABLES  COMPUTED  AT  INTERMEDIATE  POINTS  X. 
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RKGS 
• RKGS 
RKGS 
RKGS 
RKGS 
RKGS 
RKGS 
RKGS 
RKGS 
RKGS  100 
RKOS  110 
RKGS  120 
RKGS  130 
RKOS  140 
RKGS  ISO 
RKGS  160 
RKGS  170 
RKOS  180 
RKGS  190 
RKGS  200 
RKGS  210 
RKGS  220 
RKOS  230 
RKGS  240 
RKGS  250 
RKGS  260 
RKGS  270 
RKOS  280 
.RKGS  230 
RKGS  300 
RKGS  310 
RKOS  320 
RKOS  330 
RKGS  340 
RKGS  350 
RKGS  360 
RKOS  370 
RKGS  380 
RKGS  390 
RKGS  400 
RKGS  410 
RKOS  420 
RKGS  430 
RKGS  440 
RKGS  480 
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DERY  - INPUT  VECTOR  OF  ERROR  WEIGHTS . (DESTROYED) 

THE  SUM  OF  ITS  COMPONENTS  MUST  BE  EQUAL  TO  1. 
LflTERON  DERY  IS  THE  VECTOR  OF  DERIVATIVES.  WHICH 
BELONG  TO  FUNCTION  VRLUES  Y RT  A POINT  X. 

NOIM  - RN  INPUT  VALUE,  WHICH  SPECIFIES  THE  NUMBER  OF 
EQUATIONS  IN  THE  SYSTEM. 

IHLF  - AN  OUTPUT  VALUE.  WHICH  SPECIFIES  THE  NUMBER  OF 

BISECTIONS  OF  THE  INITIAL  INCREMENT.  IF  IHLF  GETS 
GREATER  THAN  10,  SUBROUTINE  RKGS  RETURNS  WITH 
ERROR  MESSAGE  IHLFsll  INTO  MRIN  PROGRAM.  ERROR 
MESSAGE  IHLF=12  OR  IHLFslS  APPEARS  IN  CASE 
PRMTC 3 3=0  OR  IN  CASE  SIGN(PRMTC'J) ) .NE .SIGN( PRMT( 2 )• 
PRMT(l))  RESPECTIVELY. 

FCT  - THE  NAME  OF  AN  EXTERNRL  SUBROUTINE  USED.  THIS 

SUBROUTINE  COMPUTES  THE  RIGHT  HAND  SIDES  DERY  OF 
THE  SYSTEM  TO  GIVEN  VALUES  X AND  Y.  ITS  PARAMETER 
LIST  MUST  BE  X.Y.OERY.  SUBROUTINE  FCT  SHOULD 
NOT  OESTROY  X AND  Y. 

OUTP  - THE  NAME  OF  AN  EXTERNAL  OUTPUT  SUBROUTINE  USED. 

ITS  PARAMETER  LIST  MUST  BE  X ,Y  .DERY . IHLF. NDIM.PRMT 
NONE  OF  THESE  PARAMETERS  (EXCEPT.  IF  NECESSARY, 
PRMT( 4 ).PRMT(S )....)  SHOULD  BE  CHANGED  BY 
SUBROUTINE  OUTP.  IF  PRMT(5)  IS  CHANGED  TO  NON-ZERO 
SUBROUTINE  RKGS  IS  TERMINATED. 

AUX  - AN  AUXILIARY  STORAGE  ARRAY  WITH  8 ROWS  AND  NDIM 
COLUMNS. 

REMARKS 

THE  PROCEDURE  TERMINATES  AND  RETURNS  TO  CALLING  PROGRAM,  IF 

(1)  MORE  THAN  10  BISECTIONS  OF  THE  INITIAL  INCREMENT  ARE 
NECESSARY  TO  GET  SATISFACTORY  ACCURACY  (ERROR  MESSAGE 
IHLFsll ) , 

(2)  INITIAL  INCREMENT  IS  EQURL  T'  OR  HAS  WRONG  SIGN 
(ERROR  MESSAGES  IHLFsi2  OR  IHLF-13). 

(33  THE  WHOLE  INTEGRATION  INTERVAL  IS  WORKED  THROUGH. 

(4)  SUBROUTINE  OUTP  HAS  CHANGED  PRMTC  5 3 TO  NON-ZERO. 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

THE  EXTERNAL  SUBROUTINES  FCT t X , Y ,DERY ) AND 

OUTPtX.Y, DERY. IHLF. NOIH. PRUT)  MUST  BE  FURNISHED  BY  THE  USER 

METHOD 

EVALUATION  IS  DONE  BY  MEANS  OF  FOURTH  ORDER  RUNGE-KUTTA 

FORMULAE  IN  THE  MODIFICATION  DUE  TO  GILL.  ACCURACY  IS 

TESTED  COMPARING  THE  RESULTS  OF  THE  PROCEDURE  WITH  SINGLE 


RKGS  460 
RKGS  470 
RKGS  480 
RKGS  490 
RKGS  600 
RKGS  510 
RKGS  520 
RKGS  630 
RKGS  640 
RKGS  660 
RKGS  660 
-RKGS  570 
RKGS  580 
RKGS  690 
RKGS  600 
RKGS  610 
RKGS  620 
RKGS  630 
RKGS  640 
.RKGS  650 
RKGS  660 
RKGS  670 
.RKGS  680 
RKGS  690 
RKGS  700 
RKGS  10 
RKGS  720 
RKGS  730 
RKGS  740 
RKGS  750 
RKGS  760 
RKGS  770 
RKGS  780 
RKGS  9Z 
RKGS  600 
RKGS  810 
RKGS  820 
RKGS  830 
RKGS  840 
■RKGS  850 
RKGS  860 
RKGS  870 
RKGS  880 
RKGS  890 
RKGS  900 
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AND  DOUBLE  INCREMENT.  RKGS  910 

SUBROUTINE  RKQS  AUTOMATICALLY  ADJUSTS  THE  INCREMENT  DURING  RKGS  920 
THE  WHOLE  COMPUTATION  BY  HALVING  OR  DOUBLING.  IF  MORE  THAN  RKGS  930 
10  BISECTIONS  OF  THE  INCREMENT  ARE  NECESSARY  TO  GET  RKGS  940 

SATISFACTORY  ACCURACY.  THE  SUBROUTINE  RETURNS  WITH  RKGS  950 

ERROR  MESSAGE  IHLF=11  INTO  MR  IN  PROGRRM.  RKGS  960 

TO  GET  FULL  FLEXIBILITY  IN  OUTPUT,  AN  OUTPUT  SUBROUTINE  RKGS  970 
MUST  BE  FURNISHED  BY  THE  USER.  RKGS  980 

FOR  REFERENCE.  SEE  RKGS  990 

RALSTON/WILF.  MATHEMATICAL  METHODS  FOR  DIGITAL  COMPUTERS.  RKGS1000 
WILEY.  NEW  YORK/LONDON.  1960.  PP.  110-120.  RKGS1010 

RKGS1020 

RKGS1030 

RKGS1040 

SUBROUTINE  RKGSIPRMT . Y ,DERY ,NDIM . IHLF .FCT .OUTP»AUX ) RKGS1050 

RKGS1060 

RKGS1070 


C 

C 

C 

C 


DIMENSION  YU)  .DERY  ( 1 ) ,RUX(  6 . 1)  .AC  4 ) .81 4 ) .C(  4 ) .PRMT(l) 

DO  1 Irl.NOIM 

RUX( 8 , 1 )s  .06666667*DERY( I ) 

X=PRMT( 1 ) 

XENDsPRMT ( 2 ) 

H-PRMK  3 ) 

PRMTC  5 )=0 . 

CALL  FCT(X.Y.OERY) 

ERROR  TEST 

IF( Hw ( XEND-X ) )38 .37 .2 

PREPARATIONS  FOR  RUNGE-KUTTA  METHOD 

A(  1 )=  .5 

A(2)=. 2928932 

A( 3 )=1 .707107 

A( 4) =.1666687 

BC  1 )=2. 

B( 2 )=1 . 

B( 3 ) = 1 . 

B(4)=2. 

C(  1 ) = .5 
C(2)=. 2928932 
CC  3 ) = 1 .707107 
C(  4 )=  .5 


C 

C 


PREPARATIONS  OF  FIRST 
DO  3 1 = 1 ,ND IM 


RUNGE-KUTTA  STEP 


RKGS1090 
RKGS1 100 
RKGS1110 
RKGS112U 
RKGS1130 
RKGS1140 
RKGS! 160 
RKGS1160 
RKGS1170 
RKGS1180 
RKGS1190 
RKGS1200 
RKGS1210 
RKGS1220 
RKGS1230 
RKGS1240 
RKGS12.S0 
RKGS 1260 
RKGS1270 
RKGS1280 
RKGS1290 
RKGS 1300 
RKGS1310 
RKGS1320 
RKGS1330 
RKGS1340 
RKGS1350 
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C 

C 

c 


c 

c 


c 

c 

c 


c 

c 

c 

c 


AUXC 1 .1 )-Y(  I ) 

RUX( 2 . 1 )=DERY(  I ) 

flUX(3.I)=0. 

RUX(6,I)-0. 

IREC=0 

HsH+H 

lHLFs-1 

ISTEPrO 

I£ND=0 


START  OF  R RUNGE-KUTTA  STEP 

4 I F ( ( XtH-XEND)*H )7 .6  »5 

5 H=XEND-X 

6 IEN0=1 

RECORDING  OF  INITIAL  VALUES  OF  THIS 

7 CALL  OUTPCX.Y.OERY.IREC.NDIM.PRMT) 
IFCPRMTC 51)40.8*40 

8 ITESTrO 

9 ISTEP=ISTEP+1 


START  OF  INNERMOST  RUNGE-KUTTfl  LOOP 
J=1 

RJ=RC J) 

BJ=B( J ) 

CJ=C( J) 

DO  11  1 = 1 ,NDIM 
R1=H*DERY( I ) 

R2=AJx( R 1-BJ*RUX(  6.1)1 

Y(I)sYU)+R2 

R2=R2+R2+R2 

RUXC  6 . 1 )=AUX( 6 . 1 ) +R2-  CJwRl 

IFC J-4J12.16.16 

J=J+1 

IFC  J— 3 513.14.13 
X=X+.6*H 

CALL  FCT(X.Y.OERY) 

OOTO  10 

END  OF  INNERMOST  RUNGE-KUTTR  LOOP 
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11 

12 

13 

14 


STEP 


IS 


TEST  OF  ACCURACY 
IFC 1 TEST  ) 16 . 16 .20 


RKGS1360 

RKGS13  0 

RKGS1380 

RKGS1390 

RKGS1400 

RkGBMlO 

RKGS1420 

RKGS1430 

RKGS1440 

RKGS14S0 

RKGS1460 

RKGS1470 

RKGS1480 

RKGS1490 

RKGS1500 

RKGS1510 

RKGS1S20 

RKGS1S30 

RKGS1S40 

RKGS1550 

RKGS1S60 

RKGS1S70 

RKGS1S80 

RKGS1590 

RKGS1600 

RKGS1610 

RKGS1620 

RKGS1630 

RKGS1640 

RKGS1650 

RKGS1660 

RKGS1670 

RKGS1680 

RKGS1690 

RKGS1700 

RKGS1  10 

RKGS1720 

RKGS1730 

RKGS1740 

RKGS17S0 

RKGS1760 

RKGS1770 

RKGS1 780 

RKGS1 790 

RKGS1800 
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I 


I'  ; 

c 

RKGS1810 

f.  i 

c 

IN  CASE  ITEST=0  THERE  IS  NO  POSSIBILITY 

FOR  TESTING  OF  ACCURACY 

RKG81820 

i. . : 

16 

00  17  1=1 .NDIM 

RKGS1830 

1 ! 

17 

RUXC4.I)=Y(I) 

RKGS1840 

i 

ITEST=1 

RKGS1860 

f ! 

ISTEP=ISTEP+ISTEP-2 

RKGS1860 

18 

IHLF=IHLF+1 

RKGS1870 

■ I 

X=X-H 

RKGS1880 

H=.5wH 

RKGS1890 

l 

DO  19  1=1 »NOIH 

RKGS1900 

» 

Y(  I )=flUX( 1 .1) 

RKGS1910 

i 

OERYC I )=flUX(2.I) 

RKGS1920 

L 

19 

flux(6.n=fiux(3.n 

RKGS1930 

f||  ; 

GOTO  9 

RKGS1940 

M 1 

c 

RKGS1960 

r-  j 

c 

IN  CASE  ITEST=1  TESTING  OF  ACCURACY  IS 

POSSIBLE 

RKGS1960 

1 ’ 

20 

IM0D=ISTEP/2 

RKGS1970 

1 

IF(  ISTEP-IMOD-IMOD321 .23.21 

RKGS1980 

1 

21 

CALL  FCTCX.Y.DERY) 

RKGS1990 

r ! 

DO  22  1 = 1 .NDIM 

RKGS2000 

I"  ; 

flUX(5.I JbY(I) 

RKGS2010 

4 

22 

AUX(7.!)=0ERY(n 

RKGS2020 

GOTO  9 

RKGS2030 

I ■ • 

c 

RKGS2040 

1 

c 

COMPUTATION  OF  TEST  VALUE  DELT 

RKGS2060 

1 

23 

0ELT=0. 

RKGS206Q 

| 

00  '24  1 = 1 .NDIM 

RKGS2070 

1 

24 

QELT=OELT+AUX( 8 . 1 ) wABSC AUXC  4 . 1 )-YC I ) ) 

RKGS2080 

IFC  DELT-PRMTC  4)128.28.25 

RKGS2090 

f >' 
; , 

c 

RKGS2100 

► 

t ' 

c 

ERROR  IS  TOO  GREAT 

RKGS2110 

; 

26 

I F ( IHLF-10126 .36.36 

RKGS2120 

26 

00  27  1 = 1 . NDIM 

RKGS2130 

■ J 

27 

AUX(4.n=Aux(s.n 

RKGS2140 

< 

IBTEP=I5TEP+ISTEP-4 

RKGS21S0 

X=X-H 

RKGS2160 

, 

IENDrO 

RKGS2170 

GOTO  18 

RKGS2180 

i'  : 

c 

RKGS2190 

c 

RESULT  VALUES  ARE  GOOD 

RKOS2200 

28 

CALL  FCT(X.Y.DERY) 

RKGS2210 

' r 

00  29  1 = 1 .NDIM 

RK0S222Q 

AUXC 1 .1  )=Y( I ) 

RKGS2230 

f 

HUXC2.I  )=DERY( I ) 

RKGS224D 

AUXC  3 . 1 )=AUX( 6.1) 

RK0S226Q 

o o o 


Y(n=oux(5.n 

RKGS2260 

29 

DERYm=fiUX(7.I) 

RKGS2270 

CALL  OUTPCX-H.Y.DERY.IHLF.NOIM.PRMT) 

RKGS2280 

IFC PRMT ( 5 ) 140.30.40 

RKGS2290 

30 

DO  31  1=1 ,NDIM 

RKGS2300 

31 

Y(n=Rux(i.n 

RKGS2310 

0ERY(I)=flUX(2.n 

RKGS2320 

IREC=IHLF 

RKGS2330 

IFC IEN0)32.32.39 

RKGS2340 

RKGS2350 

INCREMENT  OETS  D0U8LED 

RKGS2360 

32 

IHLF=IHLF-1 

RKGS2370 

l8TEP=ISTEP/2 

RK0S2380 

H=H+H 

RKGS2390 

IFC  IHLFH.33,33 

RKGS2400 

33 

IM0D=ISTEP/2 

RKOS2410 

IFC ISTEP-IMOD-IMOD )4 .34»4 

RKGS2420 

34 

IFC  0ELT-.02»PRMT(  4 1135.35.4 

RKGS2430 

35 

ihlf=ihlf-i 

RKGS2440 

ISTEP=ISTEP/2 

RKQS24S0 

H=H+H 

RKGS2460 

GOTO  4 

RKGS2470 

RKGS2480 

RKGS24&0 

RETURNS  TO  CALLING  PROORAM 

RKGS2500 

35 

iHLPall 

RKGS2S10 

CALL  FCTCX.Y.OERY) 

RKGS2520 

GOTO  39 

RKOS2530 

37 

IHLF=12 

RKGS2S40 

GOTO  39 

RKGS2550 

38 

IHLF=13 

RKGS2560 

39 

CALL  OUTPCX.Y.DERY.IHLF. NDIM, PRMT) 

RK0S2670 

40 

RETURN 

RKGS2680 

END 

RKGS2590 
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APPENDIX  C 


C . 0 ANALOG  COMPUTER  SIMULATION 

The  analog  computer  program  required  to  solve  the  three  degree-of- 
freedom  model  is  shown  in  Figure  C-l.  The  diagram  shown  is  simplified  over 
that  actually  used,  in  that  it  does  not  contain  scaling  or  gain  factor 
information,  ncr  does  it  show  certain  inverter  amplifiers  required  by  the 
non-linear  equipment.  Standard  analog  techniques  are  used  throughout  the 
program,  except  perhaps  where  multiplier  feedback  was  used  around  the  0 
and  (f>  summers  in  preference  to  resorting  to  the  use  of  less  accurate 
division  circuitry.  The  generation  and  control  of  the  firing  pulse  F(t) 
was  achieved  by  use  of  an  internal  timer,  a logic  gate  (G) , and  a delay 
flop  (P)  to  activate  two  integrators  separately  from  the  rest  of  the 
computer  equipment.  The  timer  controlled  the  firing  frequency  and  the 
delay-flop  controlled  the  duration  of  the  F(t)  function.  Equations  2-9 
required  linearization  before  they  could  be  implemented  on  the  analog 
computer.  Small  angle  approximations  were  made  and  certain  higher-order 
terms  were  removed.  The  final  set  of  differential  equations  solved  was 
as  follows: 

x = m^  L 0 + ra^.  x 02  + m^,  x <j>2  - (c^  x + k^  {x  — xQ})  - F(t) 

[I  + I + M r 2 + in  (x2  + L2)  ] 0 = m L x - (I  + m x2)  <i> 

1 cm  r cm  r r '■ 

- ir  L x f2  - c0  0 - kQ  0 - nr  g (x  - xQ)  + L F(t) 

(Ir  + mr  x2)  i - - dr  + mr  x2)  6 + \ L x 02  - - k^  +-®r  6 “ xo}  + 6 

The  main  variables  of  interest  (those  monitored  in  actual  tests)  are  out- 
put from  the  computer  solution  as  Freqqxl>  the  f°rce  felt  on  the  gunners 
shoulder,  and  the  total  pitch  angle  (0  + <S>)* 


C-l 


F(t) 


Figure  C-l  Simplified  Analog  Computer  Diagram  of 
Man-Weapon  Interaction  Equations 
C-2 
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